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INTRODUCTION 

Background 

NASA has been actively supporting the development of earth 
resources applications of passive microwave sensors over the last six 
to eight years. The driving force behind the support of passive 
microwave sensors is the high potential for estimating soil water 
interaction. Application areas that are currently receiving nwch 
attention and to which soil water is a vital requirement for their 
success are crop yield production, weather and climate modeling,^ and 
watershed management. Prior experimental work supporting the poten- 
tial of using passive microwave sensors for estimating soil water 
information has relied heavily on ground and aircraft based sensors 
[l]-[7]. Although these are natural first phase efforts, the eventual 

t 

application of such techniques will most likely be their implementa- 
tion from a space platform. 

The ground resolution cells associated with spaceborne passive 
sensors operating at low microwave frequencies are quite large because 
of limitations on the antenna size. Thus, acquisition of data from 
homogeneous, uniform areas as done with low altitude sensors will not 
be possible with spaceborne systems. At low orbit altitudes, resolu- 
tions of spaceborne passive microwave systems on the order of 5 km to 
20 km could be achieved with current technology. Resolutions of such 
dimensions will contain a mixture of the primary scene component, such 
as agricultural fields, as well as other scene components such as 
forest, urban areas, lakes, open water and rangeland. Consequently, 


results that have been demonstrated quite effectively with low alti- 
tude “high resolution" sensors might not be obtained with relatively 
coarse resolution sensors that operate from orbiting platforms. 

There are currently no passive microwave sensors In space that 
operate within the low microwave frequency range that have reasonable 
resolution dimensions. The Nimbus series of satellites have passive 
systems, but their lowest frequency is 4.9 GHz and the resolution at 
this frequency is on the order of 50 km to 60 km. The only 1.4 GHz 
passive system that has flown in space was aboard Skylab for a very 
brief period of time. Hovyever, it had a resolution of greater than 
100 km. Even with this extremely coarse resolution, encouraging 
implications were obtained when comparing the data to available soil 
water information [8]-[10]. 

Since sensor systems that have design specifications applicable 
to soil moisture estimation are not currently in orbit, it is not pos- 
sible to utilize measurements to directly determine the feasibility of 
estimating a soil water parameter from space. More importantly, if a 
system was currently in orbit, it would most likely have a fixed fre- 
quency, resolution, and look angle. It would not be possible to use 
such a systan for a thorough analysis of the effects of frequency and 
resolution on the performance of a soil moisture estimation algorithm. 

Objective 

A sensor/scene simulation program is required in order to deter- 
mine the effects of scene heterogeneity, resolution, frequency, look 
angle, and surface moisture and temperature relations on the perfor- 
mance of a spaceborne passive microwave system designed to estimate 
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soil water Information. In addition, before the expenditures for an 
orbiting passive microwave sensor system can be justified, certain 
critical questions relating to its design parameters and expected per- 
formance must be answered. The first two objectives of the project 
documented herein were to 1) develop and 2) Implement a computer pro- 
gram that could simulate the operation of a passive microwave sensor 
at L-band, C-band, and X-band frequencies for arbitrary antenna and 
orbit parameters. In addition, the model had to have the capability 
of simulating realistic scene configurations with arbitrary soil 
moisture and temperature spatial variations. A third objective was to 
utilize the model to perform soil moisture measurement feasibility 
studies. Specifically, to determine the maximum sensor resolution 
that would provide a reasonable sensitivity to soil moisture within 
the scene, and to determine the effects of scene makeup or the 
resol ut 1 on/performance relationship. 

Summary 

To meet the three objectives, the effort consisted of three basic 
tasks: 1) model determination, 2) model implementation, and 3) soil 

moisture measurement feasibility studies. Each of these tasks are 
briefly summarized below. 

Model Definition 

The primary purpose of the model is to determine the effects of 
geometric and brightness temperature variations within the scene on 
the output of a radiometer. As a result, the model consists basically 

of a digital representation of the scene with algorithms capable of 
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integrating the brightness temperature of the scene weighted by the 
one way antenna pattern. By specifying both the horizontally and 
vertically polarized brightness temperatures for each element of the 
scene, effects of the cross polarized antenna pattern can also be 
co’iputed. 

The model is defined such that the orientation of the antenna 
with respect to the scene is arbitrary. In other words, the sensor 
altitude, incident angle, and azimuth angle with respect to the scene 
are input variables. The model is also implemented such that the 
scene can be updated after each over-flight. This feature can be used 
to simulate a flight path across the scene and demonstrate the effects 
of scene geometry. It is also useful in predicting the effect of 
changing scene parameters such as soil moisture, crop type, etc., on 
the radiometer output. 

The model assumes a flat earth for both scene representation and 
flight path simula'tion. Atmospheric effects are not considered in 
this effort. However, atmospheric effects could be incorporated by 
modifying the brightness temperatures representing each element of the 
scene, or by implementing a subroutine that calculates the effect of 
the atmosphere on the brightness temperature corresponding to each 
differential solid angle involved in the antenna/scene integration. 

Model Implementation 

The model was implemented in FORTRAN for ease of implementation 
of most computer systems. The scene used in the model consists of an 
array of numbers corresponding to the percent of particular scene 
classifications contained in each pixel or array element. The array 


corresponding to the scene was generated from actual Landsat scenes 
covering a large area of south-central Texas. Each array element or 
pixel represents a 240 meter square on the earth's surface. 

The brightness temperature correspoinding to each ground scene 
pixel is calculated based upon the percentages of each scene component 
class contained therein and its relative position in the antenna 
beam. The estimated radiometer antenna temperature for each instan- 
taneous antenna resolution element within the scene is calcualated by 
integrating the weighted contribution of each pixel within the antenna 
footprint. 

Soil Moisture Measurement Feasibility Studies 

After the simulation program was tested and verified, two flight 
lines were identified across the scene and used for all simulation 
computations. These flight lines were chosen so as to provide a 
representative sampling of all possible scene configurations and com- 
ponents. Numerous simulations were computed for L-, C-, and X~band 
frequencies over these flight lines and put into a data base for 
analysis. Resolution, soil moisture, soil temperature, surface rough- 
ness, and look angle were varied between simulation runs. Also, scene 
composition in terms or the percentage of each scene class within each 
resolution element was computed and maintained by the location of each 
resolution element within the scene. 

The primary parameters investigated in this study were resolution 
or instantaneous field of view and soil moisture. The objective was 
to determine if the heterogeneity of the scene within the antenna 

resolution changes as a function of resolution size, to determine what 
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effect such a change would have on the antenna temperature soil 
moisture dependence, and to determine if there is a resolution that 
maximizes this dependence on soil moisture. In order to determine the 
effect of scene heterogeneity, soil moisture variations that would 
occur across a large scene due to precipitation patterns and soil 
composition were not addressed in this study. The effect of this 
moisture variability on the resolution requirements of an orbiting 

microwave system would be in addition to the resolution requirements 
due to scene heterogeneity. The intent was to address this question 
in the follow-on studies . 

It was originally thought that as resolution vyas decreased the 

sensitivity of the antenna temperature to soil moisture would 

increase. This was based on the assumption that the scene "purity" 

within a resolution element would improve as resolution decreased. 

This turned out to be an erroneous assumptio;.. As a result, it was 

determined that the average sensitivity of antenna temperature to 

soil moisture improves slightly as the antenna footprint size 

increased. Also, the precision (or variability) of the sensitivity 

changes as a function of resolution. Surprisingly , however, the 

highest variability occurs at middle resolutions, on the order of 20 

km. Resolutions of 5 km and 60 km have approximately equal 95% confi-^ 

dence limits on the estimate of the sensitivity to soil moisture. At 

horizontal polarization, the average sensitivity to soil moisture at 

L-band is approximately -1.5®K/percent soil moisture; at C-band 

approximately -0.85'’K/percent soil moisture; and X-band approximately 

-0.50®K/percent soil moisture. These computed sensitivities to soil 

moisture should be reasonable estimates of what can be expected from 
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an orbiting passive microwave system operating over non-mountainous 
land terrain, with the only limitation being that the pixels are less 
than 40/t forest covered. Improvements in these sensitivities could be 
expected for pixels known to consist of agricultural features. 

MODEL DEFINITION AND STRUCTURE 

In order to produce a simulation model versatile enough to be of 
use in addressing the many unanswered questions concerning the viabil- 
ity and system measurement constraints on an orbiting passive micro- 
wave system designed for soil water measurement, the system model had 
to be capable (within limits) of: 1) handling arbitrary antenna gain 
patterns, 2) integration over both like and cross polarized antenna 
gain patterns, 3) arbitrary look angle, 4) arbitrary resolution, 5) 
arbitrary altitude, 6) antenna scanning, and 7) flight path simula- 
tion. In addition, the results obtained from the simulation model are 
only as good as the scene simulated. Thus, the scene characterization 
had to be realistic in both component make up, geometry, and component 
statistics. The capability of controlling certain scene parameters 
such as roughness, soil moisture, and temperature arbitrarily over the 
scene had to be available. And, methods of estimating the brightness 
temperature of each scene pixel had to be developed as a function of 
soil moisture, soil temperature, surface roughness, microwave fre- 
quency, polarization, and emission angle. 

The simulation program generated met all of the above require- 
ments. Although all of these capabilities were not utilized in this 
study, they are expected to be needed in future efforts. Of the two 
tasks described above, 1) space platform and antenna coordinate system 
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definition/impleniGntation, and 2) ground scene definition and develop- 
ment, the latter proved to be the more difficult task and required the 
most effort. 

Space Platform Coordinate System/Inplementatlon 


The center of each ground cell within the scene is defined by its 
position in the surface-based XYZ coordinate system, Figure 1. These 
coordinates are then linked to the satellite platform coordinates 
X"'Y'"Z''* by the system of equations: 
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This transformation is accomplished by first displacing the XYZ system 
by ao, 3o* To the X’Y'Z' system. Next, the X'Y'Z' system is 
rotated about the V axis {angle 4>) to form the X''Y"Z" system. 
Finally, this latter system is rotated about the Y" axis (angle 0 ) to 
form the x"'Y"'Z"' system. The angle 4> (the first rotation of X'Y'Z' 
about Z') is defined as azimuth, while rotation about Y", angle 0 is 
defined as the roll or incident angle. 
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Antenna Temperature Equations 


The computation of the radiometer antenna temperature of that 
portion of the scene at which the antenna is looking requires an 
integration of the product of the scene brightness temperature and the 
antenna gain pattern. This integral is performed for both horizontal 
and vertical linear polarizations and is denoted by: 

RA - /Gvv(e»<{>)BTV(e.4>)dn /6vH(9.^)BTH(9.4,)dn (4a) 

JSvvTmT^ jGvH(9,^)dn 


= /GHH(Q>f)BTH(9,j»)dfl /GyH(9>({.)BTV(9,»)dQ (4b) 

^^h /GHH(9,<{))dn /Ghv( 9, ({>)dn 


where: 


BAy, BAh “ computed antenna temperature for 
vertical and horizontal polariza- 
tions, respectively. 

BTV{ 0 ,(})) - brightness temperature of the scene at 

BTH(9,(j)) vertical and horizontal polarizations, 
respectively 

Gyy(9,4>) - like polarized antenna gain patterns for 
vertical and horizontal polarizations, 
respectively 

iayH(9,(|)) - cross polarized antenna gain patterns for 

GHy{9,^) vertical and horizontal polarizations, 
respectively 

dn - differential solid angle. 


The angles 9 and 41 are identified in Figure 2. The second term in 
equations (4a) and (4b) constitute a cross polarized contribution to 
the antenna temperature. This contribution is small for good anten- 
nas. However, it is worth having the capability of computing this 
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contribution in order to determine its effect on measurements by real 
antennas that might one day be flown in space. 

Actual antenna patterns were not used in the program for the pur- 
pose of the study presented in this document. An idealized antenna 
gain pattern of the form: 


G(e,4») * G( 0) = 


sin X 
X 


f 


where 

X = tt — 

0n 

Op - half of null-to-null beamwidth of the main lobe 

For this idealized antenna pattern, there Is no cross polarized gain. 
The effect of the parameter f in the above expression is indicated in 
Table 1. 

The horizontal brightness temperatures of the scene in equations 
(4a) and (4b) are also corrected for changes in the local plane of 
incidence within the antenna footprint. If y is the angle between the 
plane determined by the sub-nadir line and the antenna axis, and the 
plane defined by the sub-nadir line and a line from the antenna and 
the pixel of interest in the scene (Figure 2), then the corrected 
scene brightness temperature is given by: 

BTV(e.4>) = TV( 0 ')cosS + TH(e')sin^Y 
BTH(e,*) = TH(6')cosS + TV(e')sinS 
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TABLE 




Effects of Parameter f on Gain Pattern Shape 


Peak Value of 
f First SI delobe 
(db down) 


3 db Beamwldth 
Compared to Null- 
to-Null Width 


0.5 

3.4 

79 

1.0 

6.7 

60 

1.5 

10.1 

50 

2.0 

13.5 

44 

2.5 

16.8 

‘10 

3.0 

20.2 

37 

3.5 

23.6 

34 

4.0 

26.9 

32 

4.5 

30.3 

30 

5.0 

33.7 

29 
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where 


I e' - local incident angle at scene pixel defined 

by 0 and ^ 

j TV(e') - vertical brightness teinperature of pixel 

^ defined by e and ^ 

I TH(6‘) - horizontal brightness temperature of pixel 

I defined by o and ^ 

j The parameters TV and TH for each type of scene pixel are computed in 

subroutines BRIGHT and CORR described below. 

I 

Integral Evaluations 

In order to evaluate the integrals in equations (4a) and (4b), 
they were transformed as follows: 


fi '■ 
i| 


il 





/G( 0 ,f)T( 0 ,^)dn 


= H T(x,y,9n,».p)G(e, 
’xy" (R‘")2 


0p dA 


(5) 


and 


fr.(»,.)rto=. y y 6(6.4.)cosendA ( 6 ) 

xy (R”')2 


where 


(R'")2 » aXaY 


and 


!i 


ii.: 

i! 


H ■ 


ir 

4 


(R'")2 - (X -«o)^ ♦ (» - 6o)^ * (Z - To)* (z) 


ir 

!!: 
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The angles 0 and <|» are given by; 


0 * COS" ^ . 11 !— (8) 

R’" 

and 

(9) 

X". 

with X'", Y*“, Z'” as defined in equations (1), (2) and (3), 

The above system of equations was implemented on the Texas A&M 
computer. The computational algorithm was designed so that for each 
nadir position of the satellite, (aj,0j), equation (4) was evalu- 
ated over the antenna footprint, Figure 3. Equations (5) and (6) were 
used to calculate the value of the indefinite integrals in equation 
(4). The footprint limit*', were established as the area included with- 
in the antenna beamwidth projected onto the ground scene. The beam- 
width of the antenna was selected as the angle included between the 
first two si delobes. By design, these two si delobes are 20 dB down 

from the gain at the center of the beam. 

The basic system algorithm accomplishes the following: 

For each X, Y, Z in the footprint, calculate 

a = COS" ^ (10) 

Ri'" 

or 

,(x 7 oi)sinosin* - (y^ - 0j)sino cos* . ^ , 

0 » cos"^ — ^ 3 II I — I— (10a) 

1 
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) 


(11) 



sin 


-1 


Y II ( 


1 


sine 


or 


‘ aj)s1n*cosG + (yi - Pj)cosocos* (Ha) 

^.= sin- Rj'sine 

where 

R,"' = [(X, - aj)2 + (y, . Bj)" + (2^ - Yj)^]V^ (12) 

The i subscript indexes each pixel in the footprint and j indexes each 
position of the satellite nadir. As these parameters are evaluated 

for each pixel in the footprint, the two sums for the indefinite 

* 

integrals, equations (5) and (6), are accumulated so that for each 
footprint, (or for each nadir position, (aj.ej)), there exists 




, JTiG,cose„Mi) 


and 


“j 


= i I ■('r“'P 


(13) 


(14) 


which are used to evaluate apparent brightness temperature for that 
footprint. 




( 15 ) 
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Similated Ground Scene 


Scene Construction 

The dominant criteria used in constructing the ground scene was 
that it must be capable of being represented digitally on magnetic 
tape or some other peripheral storage medium. Next, the scene had 
to be large enough to accommodate satellite altitudes of up to 500 km 
with antenna beamwidths of 30 degrees or less, and incident angles of 
50 degrees or less. As illustrated in Figure 4, tf?se parameters 
yield a maximum field of view of 722 km. Finally, the minimum resolu- 
tion desired was areas of approximately 40 acres. This turns out to 
be approximately 0.24 x 0.24 km; therefore, the ground scene was sized 
to be 1444 units wide. 

To insure that the ground scene was as realistic as possible, 
actual full frame classified Landsat images were used to build the 
scene. Eight Landsat images of central and east Texas classified into 
various vegetation, water, and urban classes by the Texas Parks and 
Wildlife Department were used. Figure 5 illustrates the area covered 
by the eight Landsat scenes. Appendix A contains the classified 
images as produced by the Texas Parks and Wildlife Department, along 
with descriptions of what the classifications mean. 

Each classified Landsat image contained a maximum of 1824 by 2048 
Landsat pixels. Further, each classified map had different classes. 
There was inconsistency in the definition format of the scene classes 
on different maps (i.e., the digital count representing grass would 
appear as a different value on each tape containing the classified 
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FIGURE 4. Limitations of the simulation model in terms of field of view 
requirements and scene size. 
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7 BRYAN 

8 HOUSTON 



FIGURE 5. 


Areas covered by the eight Landsat scenes used to generate the 
simulated scene. 





data from each scene). In order to use the Texas Parks and Wildlife 
classified Landsat images to produce a ground scene, the 80 m by 80 m 
Landsat pixels had to be aggregated to decrease the number of scene 
pixels. Also, the classes as defined by the Texas Parks and Wildlife ■ 
Department had to be reassigned into classes meaningful to microwave 
emission computations. The class assignments in the Texas Parks and 
Wildlife classification were analyzed and transformed into the follow- 
ing six classes described in Table 2. The procedure for performing 
this reduction reclassification is described below. Models for com- 
puting the brightness temperature of each of the classes defined in 
Table 2 as a function of soil moisture, temperature, and roughness are 
also described below. 

Converting Classified Landsat Types to TAMU Percent Class Types 

Each pixel in the simulated scene was created by aggregating 3 by 
3, 80 m by 80 m pixels in the classified Landsat maps to one 240 m by 
240 m pixel in the simulated scene. In order to maintain the class 

composition of each aggregated 240 m by 240 m pixel, the percentage 
coverage of each class within the aggregated pixel was computed and 
maintained. These percentages were used to compute the proper bright- 
ness temperature for the mixture of classes within the aggregated 
pixel . 

The Texas Parks and Wildlife tapes (MAPTAP) were formatted as 
scan lines of up to 2048 pixels per line such that each byte (8 bits) 
represented an 80 x 80 meter ground cell. Up to 1824 lines made up a 
single scene. The absolute binary value of each byte on the scan line 
was indicative of a land cover class (as defined by the Texas Parks 


Class 

Description 

1 

Water 

2 

Bare Soil 

3 

Urban 

4 

Mixed Soil and vegetation 

5 

Fully vegetated 

6 

Forest 


and Wildlife Service) within the respective ground cell. For example, 
in the Kerrville scene, MAPTAP No. 0068, the byte va jes range from 0 
to 24 (decimal) as indicated in the legend for that scene in Appendix 
A. To convert these tapes to a more usable size and format, a 
reduction/reclassification transformation was performed. As indicated 
in Figure 6, each 6 x 6 set of bytes each byte representing one 
Landsat pixel on the MAPTAP was reduced into a set of 2 x 2 words (16 
bits) of percent class data. Each word represents one 240 m x 240 m 
scene pixel. 

Each MAPTAP was scanned to determine the percentage of the six 
pre-selected classes identified above that existed within each 3x3 
Landsat pixel set. Since the MAPTAP vegetation classes did not match 
the 6 classes shown above and at the bottom of Figure 6, a selective 
transformation was made in the MAPTAP data. For example, on the 
Kerrville scene all MAPTAP values of 1, 2, 3, or 4 were considered to 
be TAMU Class 5; i.e. , Fully Vegetated. The selected transformation 
that was used for each of the eight scenes varied according to the 
MAPTAP value assignments. Since each tape was different, a separate 
transformation scheme was selected for each tape (scene). The 
right-hand column of Appendix A shows the groupings used for each of 
the Texas Park and Wildlife scenes used. 

Each 240 m by 240 m simulated scene pixel in the TAMU class data 
tape was represented by one 16 word bit. Each octal word within the 
16 bits contained the percent of class coverage within that aggregated 
pixel for one of the six classes. In this manner, the 16 bit word 
contained the percentage coverages for five of the six classes. The 
percentage of the sixth class (forest) was implied by subtracting the 
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3 words (6 bytes) 





\ 
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1 

1 

1 



1 

1 

1 




1 

1 



1 

1 

1 



1 

1 

1 



1 

, 

1 








Texas Parks and Wildlife MAPTAP 
format. Each byte 255 counts) 
represents a Land cover class and 
corresponds to a Landsat 80 x 80 
meter resolution element 



I 


Each 6 X 6 bytes of MAPTAP infor- 





mation was mapped into 2x2 


1 

1 

words. Each word represents per- 


1 

1 

cent class data for a 240 x 240 




meter ground cell. Explanation of 


i \ how each octal word represents 

2 words * class assignments is in Figure 7, 


.T ' 

I 
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CLASS ASSIGNMENTS IN THE “CLASS DATA" TAPE 


1 - 

WATER 

4 - MIXED SOIL & VEGETATION 

2 - 

BARE SOIL 

5 - FULLY VEGETATED 

3 - 

URBAN 

6 - FOREST 


FIGURE 6. Procedure for aggregating and mapping the Texas Parks and 
Wildlife Landsat classifications into the class used in the 
simulation program. 
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sum of the other five from 100. Encoding the percent class data into 
the 16 bit word describing each pixel of the simulated ground scene is 
illustrated in Figure 7. The [ means the integer equivalent 

of the floating point expression. 

Another problem was that the map tapes made no distinction 
between bare soil and urban area. Since these two areas have dis- 
tinctly different brightness temperatures, adjustments had to be 
made. First, known urban areas were identified manually and their 
coordinates were given to the computer. These areas were then auto- 
matically changed from bare soil to urban as they were encountered 
during the transformation. Next, to get a more general mix of urban 
over the entire scene, whenever a 3 by 3 pixel set in the vegetation 
tapes was found to be fully bare soil, it was re-defined as urban. 
After each transformation, the resulting class map was plotted in 
grey-scale format to verify that an appropriate distribution had been 
achieved. 

Using corner coordinates supplied with the Texas Parks and 
Wildlife sub-scene map tapes, the entire class data ground scene was 
grouped into one tape file as illustrated in Figure 8. The tape files 
contain 1650 records of 2496 words each. The process of constructing 
the composite class data scene was a nine step sequence. First, the 
entire scene was filled with Class 6 (forest) as a background. Each 
sub-scene of class data was over-layed starting in the upper right 
hand corner with the Lufkin sub-scene, and continuing with Houston, 
Corsicana, Bryan, Waco, Austin, Brownwood, and finally Kerrville. The 
entire class scene was then grey-scaled for each class. These maps 
are included in Appendix B. 


25 



ORIG'NM- 
Of POOR 


16 Bit word corresponding to one 240 x 240 m simulated scene pixel 


0 1-3 


4-6 


7 - 9 


10-12 13 - 15 


CLASS 5 


- - 

' (f 7 
L9. /.Ji 


*4096 


CLASS 4 


■ (f^)* i 

L^* J I 


CLASS 3 


— 

= (|^)* 7. , 

L®- J I 


Notes: 


1. Ei is the number of 
Landsat pixels in the 
3x3 array of Landsat 
pixels corresponding to 
the ith class, 

2. The subscript I indicates 
that the number in brackets 
is an integer. 


CLASS 2 -■ 


= fe)* 7] , 

L- * J I 


CLASS 1 


= Hr)* 7. 

L J I 


CLASS 6 IMPLIED 


FIGURE 7. Construction of the 16 bit word that describes the percentages 
of classes in each 240 x 240 meter simulated scene pixel. 
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Scene Utilization 


Use of the composite ground scene and the related computer 
software has specific limits as implied in Figure 9. The altitude 
limits are related to antenna beamwidth. The scene pixels represent 
approximately 14 acres each. The entire scene represents a surface 
area of 396 by 599 km. 

The system initialization and operation concept is illustrated 
in Figure 10. It is important to understand that the starting point 
is defined by the initial coordinates of the antenna beam footprint 
(center line); i.e., (Xq.Yo)* This, along with altitude, incident and 
azimuth angles, define the initial system geometry. Also required are 
beamwidth and displacement rates. During this phase of the project, 
zero rates were used for all variables except the nadir coordinates. 
Satellite movement was accomplished by stepping the nadir coordinate 
from (ao*3o) To (07,8^) by the rates Aa, Ap. 

Since the initial position of the satellite is specified in 
terms of the antenna beam center line on the surface, the initial 
nadir point must be evaluated using the concept illustrated in Figure 
11. At the beam center line, 

X"' = Y'" =0 (16) 

and 

Z"’ =-^ (17) 

cos« 

then 

Rj = Z"'sin 0 = -Htano (18) 

so that 
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FIGURE 9. Major simulation parameters that must be specified properly in 
order to maintain consistency between the antenna footprint 
size and the simulated scene size. 


I 
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FIGURE 11. Calculation of the initial nadir point. 
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Rssin* = ao - Xq 

(19) 

Rscos* = Yq - Bo 

(20) 

®o ~ RsSin*+ Xq 

(21) 

Bo Vq - Rscos* 

(22) 


Now, knowing the starting condition it remains to determine 
reasonable limits on the integration over the antenna footprint. 
These limits were selected as upper and lower boundaries on a 
rectangle which always contains the smaller rectangle whose sides are 
parallel to the major and minor axis and tangent to the sides of the 
footprint ellipse, Figure 12. The lower and upper limits on X and Y 
are derived as follows (Figure 12); 


a = -H tano - c (23) 
b = (H/coso)tan 8W/2 (24) 
c = H tan (o “ BW/2) (25) 
d = H tan (o + BW/2) - a - c (26) 

then 

YU = + d»cos 4 > + b.sin* (27) 
YL = Yi + a. cost - b.sin* (28) 
XU = Xj + a. sin* + b.cos* (29) 
XL = Xi + d»sin* - b.cost (30) 


Subsequent program steps in calculating the brightness tempera- 
ture over the antenna footprint are outlined in the Computational Flow 
Chart in Appendix C. For each nadir position defined by (ao + kAa, Bo 

32 


origima'- page ra 

OF POOR QUAUTT 



FIGURE 12. Definition of the parameters describing the limits of the 
antenna footprint integration. 
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+ an apparent brightness temperature is calculated. Thus, the 

frequency of the data points calculated, as the satellite passes over 
the scene, and the direction of the movement across the scene, are 
governed by the selection of ao and ap. The index k simply increments 
from zero until the edge of the footprint encounters an edge of the 
composite scene. 

The software system was implemented using the concept outlined 
in Figure 13, The main program was run either on the TI980 mini- 
computer or on the AHDAHL470 (IBM 360 Operating System). Coding was 
all in FORTRAN. Initialization was through the card reader or key- 
board (T1980). Class data (ground scene) was on magnetic tape. The 
brightness temperatures and other associated data were output both in 
tabular form and on magnetic tape (compressed data set). All plots 
were made on the T1980 with output to a VERSATEC matrix printer. 

All software code is included in Appendix D. These include the 
following: 

STEP 1 PROGRAM - Maps vegetation 'type' k ata for each of 8 sub- 
scenes to each of 8 subscenes of '% class' data 

sets. Each unique mapping algorithm is 

included.’ 

FILL DATA - Fills composite scene with forest background 
prior to performing individual overlays. 

STEP 2 PROGRAM - Performs overlay operation; i.e., placing each 

subscene onto the composite scene at the 

correct position. 

SATELLITE MODEL - Accepts initialization data and performs 

simulation of moving satellite over the scene, 
producing brightness temperature estimates for 
each position of the nadir point. 
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COMPUTATION 


The models used to compute the brightness temperature of each 
pixel in the simulated scene (BTH and BTV in equation (4)) were devel- 
oped under subcontract to the University of Houston at Clear Lake 
City. 

Objective and Approach 

The objective of the research performed under the subcontract 
was to develop a mathematical model that can be used on a computer 
system to simulate passive microwave radiometer measurements over 
heterogeneous earth scenes under varying soil moisture, cultural fea- 
ture, land cover, and atmospheric conditions. The algorithms that 
make up the mathematical model are limited to operating frequencies of 
1.42 GHz (L-band), 4.8 GHz (C-band), and 10.7 GHz (X-band), and gen- 
eral land cover types of water, forest, grassland, partially vege- 
tated, bare soil, and urban classes. The mathematical models predict 
the microwave emission of land or water surfaces for each of the six 
basic land cover types (i.e., water, fully vegetated (grassland), 
forest, mixed soil and vegetation, bare soil, and urban as specified 
in Table 2). The predicted emissions were designed to be valid at a 
nadir viewing angle of 50 degrees and for both horizontal and vertical 
polarization. Microwave emission at viewing angles other than 50 
degrees was obtained through the use of yet another mathematical 
model . 

To make the various computer models as realistic as possible, 
the class emission models were based upon empirical measurements with 
truck -mounted or airborne microwave radiometers where possible. In 
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most cases, such measurements existed and were documented in the open 
literature. Some gaps were evident, however. In these cases, theore- 
tical models were used. As a matter of procedure, the information and 
models developed at the University of Houston at Clear Lake City were 
sent to the Remote Sensing Center throughout the contract period. The 
actual encoding of the algorithms on a computer system was done at the 
Remote Sensing Center at Texas A&M University. 

The various mathematical models are described in detail in the 
following sections. 


Subroutine BRIGHT 

To provide a realistic data base of land cover types for the 
simulation program, as previously described, investigators at the 
Remote Sensing Center obtained a set of land cover themes based on 
classified Landsat images for the general area of central and eastern 
Texas. The area had been classified into several land cover catego- 
ries which in turn could be assigned to one of the following six 
classes: open water, vegetation (total cover), bare soil, mixed vege- 
tation and bare soil, forest, and urban. These data were given in a 
data base where each data cell represented one Landsat picture element 
having a spatial area of about 1.1 acre. It was determined by members 
of this current research effort that the given data base was too large 
to handle. So, it was decided to reduce the data base by combining 56 
Landsat picture elements into one composite area. This procedure 
results in the basic information unit for land cover being heterogene- 
ous, that is, many classes can be represented in the 40 acre data 


cell . 
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It was necessary to preserve some aspect of the mixture of 
classes within the composite data cell. Since a 16 bit word was 
available to encode the class type, it was decided to include some 
information about the percentage area of each class in the single data 
word used to represent the composite area in the reduced data base. 
After the leading sign bit, the data word was made up of five three- 
bit word slices. Each three-bit slice can be interpreted as follows: 
Bit Pattern Percent of Area Covered 


000 

001 

010 

Oil 

100 

101 

no 

111 


0.0?J 

18.75 

31.25 

43.75 

56.25 

68.75 

81.25 

100.00 


For example, suppose that 255i of the picture elements represented 
bare soil, 15% represented water, and the rest (60%) was mixed 
vegetation and soil. From left to right, the slices represent 
vegetation, mixed soil and vegetation, bare soil, urban, and open 
water. The class, forest, is included as a residual. The above 

percentages would then be encoded as follows (after the sign bit): 
000100001000001. In turn, this particular pattern would be 
interpreted as follows: 


Vegetation 0.0% 

Mixed Soil & Vegetation 56.26 

Bare Soil 18.75 

Urban 0.0 

Open Water 18.75 

Subtotal 93.75 


A residual of 6.25% would be assigned to the class of forest 
in this example. The example clearly shows the fact that errors exist 
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in the categorical way of expressing the percentage mixture of each 
subclass. The advantage gained is the increase in the speed of 
computation with a data base that has been reduced by a factor of 40. 


FORTRAN Code for Subroutine BRIGHT 


The FORTRAN code for Subroutine BRIGHT is given in Appendix 
0. The purpose of Subroutine BRIGHT is to calculate the horizontally 
and vertically po'^arized brightness temperature for a given composite 
data area at a viewing angle of 50 degrees from the nadir for given 
land cover types and conditions. The call for Subroutine BRIGHT is as 
fol lows: 

CALL BRI6HT(SM,ICW,NB,TP,ROU,BTV50,BTH50} 

The elements of the call are defined as follows: 

SM - Soil moisture parameter taken to be the volumetric soil 
moisture content in percent for Miller Clay. A field 
capacity of 100% is reached when SM==38. 

ICW - A 16-bit integer binary word that contains information 
about the distribution of ground cover types within the 
composite area cell, (see above for code) 

NB - Frequency band number (NB=1 for L-band, NB=2 for 
C-band, and NB=3 for X-band) 

TP - Temperature parameter (the temperature that dry, bare 
soil would have under the same weather conditions, 
given in degrees Celsius) 

ROU - Roughness parameter (used for bare soil and mixed 
vegetation and soil calculations -- use ROU=0.0 for 
smooth fields and R0U=0.6 for rough soils) 

BTV50 - Brightness temperature of the composite area cell for 
vertical polarization and for a zenith angle of 50 
degrees (degrees Kelvin) 

BTH50 - Brightness temperature of the composite area cell for 
horizontal polarization and for a zenith angle of 50 
degrees (degrees Kelvin) 
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Subroutine BRIGHT uses several other subroutines in the pro- 
cess of the calculation of the brightness temperatures of the compo- 
site cell. Also, the information encoded in the parameter, ICW, is 
is unpacked and used in the calculation. Subroutine BRIGHT uses the 
following subroutines: 


Subroutine WATER 
Subroutine URBAN 
Subroutine BARE 
Subroutine MIX 

Subroutine VEG 
Subroutine FOREST 
The outputs from Subroutine 
used in the Subroutine BCORR which 
temperature components at angles 
zenith. 

The other subroutines menti 
in the following sections. 


Emission of open water 

Emission of urban areas 

Emission of bare fields 

Emission of mixed vegetation and 
soil 

Emission of grassland 
Emission of forest 
BRIGHT (i.e., 8TV50 and BTH50) are 
estimates the polarized brightness 
other than 50 degrees from the 

oned above are discussed in detail 


Brightness Temperature as Function of Angle 
(Subroutine BCORR) 


Based upon the many measurements of the brightness temperature 
of various objects as reported throughout the literature, a simple 
extrapolation scheme was adopted to extend calculated brightness 
temperatures to angles of viewing needed in the program calculations. 
A program was developed to calculate the emissivity of a smooth, 
homogeneous dielectric material interfacing air. The calculations 
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were perfoniied for viewing angles (measured from the normal to the 
interface surface) from 0 to 80 degrees in steps of one degree. The 
calculations were performed for a number of material types (liquid 
water, dry soil, wet soil). It was noted that the normalized 
distribution of emissivities with viewing angle was little affected by 
the choice of dielectric constant. Thus, it is possible to model the 
angular distribution with a tabular function as follows: 

BTH = BTH50 + FH*0ELBT 
BTV = BTH50 + FV*DELBT 


where 

BTH “ brightness temperature at some given angle for horizontal 
polarization 

BTV - brightness temperature at some given angle for vertical 
polarization 

FH “ form factor for horizontal polarization which is a 
function of viewing angle (table lookup) 

FV - form factor for vertical polarization which is a function 
of viewing angle (table lookup), 

DELBT « difference between the brightness temperatures (vertical 
and horizontal polarization) at a viewing angle of 50 
degrees, i.e., 

DELBT = BTV50 - BTH50 

BTH50 - calculated brightness temperature for horizontal 

polarization and for a viewing angle of 50 degrees as 

given from Subroutine BRIGHT. 

BTV50 - calculated brightness temperature for vertical 

polarization and for a viewing angle of 50 degrees as 

given from subroutine BRIGHT 

The values for FH and FV adopted in the present algorithm are 
given in Table 3 and plotted in Figure 14. 

The appropriate value for FH and FV is obtained by conversion of 
the angle to an integer index value. Using the index value, the 
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TABLE 3, Values for the form factors, FH and FV, as a 
function of nadir angle. 


Angle (degrees) 

FH 

FV 


0 

0.540 

0.540 


1 

0.540 

0.540 


2 

0.539 

0.540 


3 

0.538 

0.541 


4 

0.537 

0.542 


5 

0.535 

0.544 


6 

0.534 

0.546 


7 

0.531 

0.548 


8 

0.529 

0.551 


9 

0.526 

0.554 


10 

0.522 

0.557 


11 

0.519 

0.561 


12 

0.515 

0.565 


13 

0.510 

0.569 


14 

0.506 

0.574 


15 

0.500 

0.579 


16 

0.495 

0.584 


17 

0.489 

0.590 


18 

0.483 

0.596 


19 

0.476 

0.602 


20 

0.469 

0.609 


21 

0.461 

0.617 


22 

0.453 

0.624 


23 

0.445 

0.632 


24 

0.436 

0.641 


25 

0.426 

0.649 


26 

0.417 

0.659 


27 

0.406 

0.668 


28 

0.395 

0.378 


29 

0.384 

0.689 


30 

0.372 

0.699 


31 

0.360 

0.711 


32 

0.347 

0.722 


33 

0.333 

0.734 


34 

0.319 

0.747 


35 

0.305 

0.759 


36 

0.289 

0.773 


37 

0.274 

0.787 


38 

0.257 

0.801 


39 

0.239 

0.815 


40 

0.222 

0.830 
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1 ^ 

TABLE 3. (con't). Values 

for the form factors, 

FH and FV, as 



a tabular function of nadir angle. 


1. 

Angle (degrees) 

FH 

FV 



41 

0.203 

0.845 



42 

0.184 

0.861 


r 

43 

0.164 

0.877 


1 __ 

44 

0.143 

0.894 



45 

0.121 

0.911 


1 - 

46 

0.099 

0.928 


' 

47 

0.075 

0.945 



48 

0.051 

0.963 



49 

0.026 

0.982 



50 

0.000 

1.000 


51 

-0.027 

1.019 



52 

-0.055 

1.038 



S3 

-0.084 

1.057 


■t, 

54 

-0.114 

1.076 



55 

-0.145 

1.095 


! 

56 

-0.177 

1.115 



57 

-0.211 

1.134 



68 

-0.245 

1.153 



59 

-0.281 

1.172 



60 

-0.3J^ 

1.190 



61 

-0.356 

1.208 



62 

-0.396 

1.225 


g 

63 

-0.437 

1.241 



64 

-0.480 

1.256 



65 

-0.523 

1.270 


2 

66 

-0.569 

1.282 


i 

» 

67 

-0.616 

1.292 



68 

-0.664 

1.300 



69 

-0.714 

1.305 


k 

70 

-0.766 

1.306 


^ - 

71 

-0.820 

1.304 



72 

-0.875 

1.297 


g 

f 

73 

-0.932 

1.285 


f. 

t •: 

74 

-0.991 

1.266 



75 

-1 .051 

1.239 


r ■ ' 

76 

-1.114 

1.203 


\ 

77 

-1.179 

1.157 



78 

-1.246 

1.098 



79 

-1.314 

1.025 


1 

80 

-1.385 

0.933 


p 

H 

n 
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FIGURE 14. Plot of the form factor used to estimate brightness temperature 
as a function of incident angle. 
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needed value for FH and FV is looked up in the table of values of FH 
and FV. The brightness temperatures computed in subroutine BRIGHT at 
50® (BTV50 and BTH50) are computed to the appropriate incident angle 
in subroutine BCORR. The listing of this subroutine is contained in 
Appendix D. The call for subroutine BCORR is as follows: 

CALL BCORR ( BTV50 , BTH50 , TN . BT V , BTH ) 

The elements of the call are determined as follows: 

TN - incident angle at the pixel of interest 
BTV - vertical brightness temperature at angle TN 
BTH - horizontal brightness temperature at angle TN 

Subroutine WATER 

The emissivity of water at microwave frequencies is a function of 
water temperature, water salinity, frequency, angle of viewing, and 
polarization Linear relationships were derived between emissiv- 

ity and water temperature at a 50 degree incident angle for both hori- 
zontal and vertical polarization based upon calculations made by Paris 
[2] in which the water surface was modeled as being a flat, homgeneous 
dielectric material overlain by air. These relationships were 
derived for the three frequency bands of interest to this research 
(i.e., L-, C-, and X-band). The results of the calculations are as 
follows: 

L-band: 

EH50 = 0.256 + TWC*0.000467 
EV50 = 0.505 + TWC*0. 000767 

Note: EH50 is the horizontal emissivity at 50 degrees 

viewing angle, EV50 is the vertical emissivity at 
50 degrees viewing angle, and TWC is the water tem- 
perature in degrees Celsius. 
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C-band: 


EH50 = 0.265 
EV50 = 0.522 


X-band; 

EH50 = 0.288 TWC*0.0003 

EV50 = 0.557 - TWC*0.0005 

With the value of the emissivity, the brightness teinperaure is calcu- 
lated by the following formulas: 

BTVW = EV50*TWK + (1. - EV50)*SKYT 

8THW = EH50*TWK + (1. - EH50)*SKYT 

where 

TWK - temperature in degrees Kelvin 

SKYT - sky brightness temperature at 50 degrees 

viewing angle (SKYT = 6.0 for L-band, 8.0 for C-band and 
10.0 for X-band) 

BTVW - value of the brightness temperature of the water for ver- 
tical polarization and for a viewing angle of 50 degrees, 

BTHW - the value of the brightness temperature of the water for 
horizontal pol ari z.ition and for a viewing angle of 50 
degrees. 


The FORTRAN code for Subroutine WATER is given in Appendix D. 

The call for Subroutine WATER is as follows: 

CALL WATER (NB,TP,BTN,BTHW) 

The elements of the call are defined as follows: 

NB - frequency band number (as before) 

TP - temperature parameter (as before) 

BTVW - brightness temperature of water for veritcal polarization 
and a viewing angle of 50 degrees 
BTHW - brightness temperature of water for horizontal 
polarization and a viewing angle of 50 degrees. 
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Subroutine URBAN 


An ‘Urban area is an extremely complex environment for which 

calculation of the brightness temperature is required. Water can pond 
on the tops of buildings. Many different types of materials are found 
in such areas. Also, natural vegetation and forest are mixed in such 
areas. A very simple approach was adopted in the face of such 

complexity. It was assumed that the emissivity of the urban area over 
all was simply the following: 

EV50 = 0.96 
EH50 = 0.86 

Moreover, it was assumed that the temperature of the urban scene 
was the same as that of dry, bare soil under the same climatic 

conditions. 

The assumed values for the emissivity of an urban area were 
based upon measurements over asphalt by Straiton and Talbert [13], 

over asphalt, gravel road, and other manmade materials by Straiton ^ 
al . [14], and on other such materials by Porter and Florance [15]. 

The FORTRAN code for Subroutine URBAN is given in Appendix D. 
The call for Subroutine URBAN is as follows: 

CALL URBAN(TP,BTVU,BTHU) 

The elements of the call are defined as follows: 

TP - temperature parameter (as before) 

BTVU - brightness temperature of an urban area for vertical 
polarization and for a viewing angle of 50 degrees 

BTHU - brightness temperature of an urban area for horizontal 
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Subroutine BARE 


The microwave emission of bare soil has been under study for 
many years by investigators in the Joint Soil Moisture Experiment 
sponsored by NASA. In general, they have found that the emissivity 
of soil is a function of soil moisture in the upper layers, soil 
roughness, soil type, angle of viewing, polarization, and frequency. 
The brightness temperature of soils is influenced by soil temperature 
as well as the forementioned parameters. 

For the purpose of the algorithm reported in this report, the 
time of day is taken to be about 2 p.m. local time, a time when the 
surface of the soil is near its maximum temperature. For a given soil 
moisture condition, the microwave radiation that emerges from the soil 
surface is emitted from a zone of soil bounded on the top by the air- 
soil surface. The thickness of that zone is determined by soil mois- 
ture content itself and the frequency of the radiometer. Also, the 
relationship between emissivity and soil moisture is non-linear. In 
this algorithm, the non-linear function is approximated by two linear 
functions that are continuous. One function covers the range of soil 
moisture from 0 to 12% (by volume), and another function covers the 
range of soil moisture greater than 12% (by volume). 

The most complete study of the microwave emissive properties of 
a soil type is that of Miller Clay near the Remote Sensing Center, 
Texas A&M University, College Station, Texas. The main results of 
several summers of measurements from truck-based microwave radiometers 
have been reported by Newton [4]. Recently, Choudhury el^ [11] 
found that the roughness effects on a field of bare soil can be 
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modeled by an exponential function* These results have been incorpor- 
ated into the algorithms below. 

Soil moisture itself has an effect on soil temperature. This 
fact is included in the algorithms proposed in this study. Figure 15 
demonstrates the manner in which this is taken into account in the 
simulation for a temperature parameter (TP defined below) of 40®C. 


X-band Algorithm 

At X-band, the emission occurs from the uppermost layer for all 
ranges of soil moisture. The effect of soil moisture on soil 

temperature (T6K in degrees Kelvin) is given as: 

TGK = 273.15 + TP - 0.87*SM 


where TP is the temperature parameter (degrees Celsius), and SM is the 
soil moisture (percent by volume). For Miller Clay, SM ranges from 
zero to 38% at field capacity. Moisture contents greater than field 
capacity are possible. If SM is greater than 38%, TGK is set equal to 
TP plus 240.15. 

The equations for the emissivity are as follows: 


SM < 12%: 

EH50 = 0.91 - 0.00917*SM 
EV50 = 0.99 - 0.0025*SM 
SM > 12%: 

EH50 = 0.96 - 0.0135*SM 
EV50 = 1.05 - 0.0077*SM 
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C-band Algorithm 


Little work has been done at C-band In passive microwave remote 
sensing of soils. It is assumed here that the effects of soil mois- 
ture lie halfway between the effects at L-band and at X-band. The 
main departure is the modeling of dielectric constant from which emis- 
sivities are gathered. The algorithm at C-band is as follows: 

TGK = 260.15 - 0.53'^SM + TP 
SM < 12%: 

EH5U = 0.86 - 0.00833'^SM 
EV50 = 0.97 - 0.0025*SM 
SM > 12%: 

EH50 = 0.92 - 0.0135*SM 
EV50 = 1.04 - 0.00846*SM 

L-band Algorithm 

The algorithm at L-band is as follows: 

TGK = 250.15 + TP - 0.26*SM 
SM < 12%: 

EH50 = 0.90 - 0.00917*SM 
EV50 = 0.98 - 0.0025*SM 
SM > 12%: 

EH50 = 0.96 - 0.0139*SM 
EV50 = 1.047 - 0.00808*SM 

FORTRAN Code for Subroutine BARE 

The FORTRAN code for Subroutine BARE is given in Appendix D. 
The call for Subroutine BARE is as follows: 
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CALL BARE(NB,TP,SM,ROU,BTVB,BTHB,RH,RV) 


The elements of the call are defined as follows: 


NB - Frequency band number (as before) 

TP - Temperature parameter (as before) 

SM - Soil moisture parameter -- the soil moisture 

for a Miller Clay soil in percent by volume. For other 
soils, translate the soil moisture to that of Miller Clay 
by use of wilting point (12^) and field capacity (38%). 

ROD - Roughness parameter (as before) 

BTV8 - Brightness temperature for bare soil for vertical polari- 
zation and for a viewing angle of 50 degrees. 

BTHB - Brightness temperature for bare soil for horizontal 
polarization and for a viewing angle of 50 degrees. 

RH - Reflectivity for horizontal polarization and for a view- 
ing angle of 50 degrees (parameter is needed for Subrou- 
tines MIX and VEG) 

RV - Reflectivity for vertical polarization and for a viewing 
angle of 50 degrees 


Roughness Algorithm 


A rough surface has a higher emissivity than does a smooth 
surface. The algorithm used to account for the effect of surface 
roughness is as follows: 

RH = 1.0 - EH50 
RV =1.0 - EV50 
RFAC = EXP (-R0U*0 .4132) 

RHnew = RHold*RFAC 
RVnew = RVold*RFAC 
EH50new = 1.0 - RHnew 
EVBOnew = 1.0 - RVnew 


52 


Then, to complete the calculation for BTHB and BTVB: 

6THB = EHBO^TGK 
BTVB = EH50*TGK 

An example of this effect is seen in Figure 16. 

Subroutine VEG 

In this subrputine, it is assumed that a 100 percent grass type 
vegetal cover exists on the land. In such cases, the temperature of 
the vegetal canopy is moderated by the processes of evaporation and 
transpiration* The algorithm below attempts to allow for this 

effect. Figure 17 illustrates the dependence of the plant canopy 

I 

temperature on the temperature parameter. 

The brightness temperature of vegetation is computed using the 
reflectivities and brightness temperatures obtained from Subroutine 
BARE under smooth surface conditions (i.e., R0l)=0.0). The 
reflectivity is modified according to a computed vegetation factor, 
VFAC, which depends upon the soil moisture content. This factor is 
based upon measurements made by Newton [4]. The FVAC reduces the 

computed reflectivity. The vegetation factor is applied fully for 

L-band calculations, partially for C-band calculations, and is not 
used for X-band calculations. This means that at X-band, no 

penetration of the vegetal canopy occurs. In the case of X-band, a 
constant set of emissivities is assumed, viz., EH50=0,92 and 

EV50=0,95. The L-band algorithm is given below: 

TVC - (TP - 25.0)*0,24 + 25.0 (Correction for canopy 

temperature) 

TVK = TVC + 273.15 

CALL BARE(NB, TVC, SM,0.0, BTVB, BTHB,RH,RV) 
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FIGURE 17. Example of how the ten^erature parameter is assumed to effect 
the water or plant canopy temperature. 


VFAC = 0.8 - 0.00395*SM 


XRH = RH^VFAC 

XRV = RV*VFAC 

The FORTRAN code for Subroutine VEG is given in Appendix 0. The 
call for Subroutine VEG is as follows: 

CALL VEG(NB,TP,SM.8TVV,8THV) 

The elements of the call are defined as follows: 

NB - Frequency band number (as before) 

TP “ Temperature parameter (as before) 

SM - Soil moisture parameter (as before) 

BTVV - Brightness temperature of the vegetal canopy for vertical 
polarization and for a viewing angle of 50 degrees 

BTHV - Brightness temperature of the vegetal canopy for horizon- 
tal polarization and for a viewing angle of 50 degrees. 

Subroutine NIX 

Subroutine MIX is a subroutine that calculates the brightness 
temperatures for a partially vegetated field. Subroutines BARE and 
VEG are called in the program. The results are averaged together to 
get the mixed result. 

The FORTRAN code for Subroutine MIX is given in Appendix D. The 
call for Subroutine Mix is a follows: 

CALL MIX(NB,TP.SM,ROU,BTVM,BTHM) 

The elements of the call are defined as follows: 

NB - Frequency band number (as before) 

TP - Temperature parameter (as before) 

SM - Soil moisture parameter (as before) 

ROU - Roughness parameter (as before) 
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BTVM - Brightness temperature for the partially vegetated field 
for horizontal polarization and for a viewing angle of 50 
degrees 

BTHM - Brightness temperature for the partially vegetated field 
for horizontal polarization and for a viewing angle of 50 
degrees. 

Subroutine FOREST 

Subroutine FOREST simply treats the surface as a thick vegetated 
canopy. Subroutine VEG is used with NB-3 regardless of the band used. 

The FORTRAN code for Subroutine FOREST is given in Appendix D. 
The call for Subroutine FOREST is as follows; 

CALL FOREST(TP.BTVF,BTHF) 

The elements of the call are as follows: 

TP - Temperature parameter (as before) 

BTVF - Brightness temperature of the forest for vertical polari- 
zation and for a viewing angle of -SO degrees 

BTHF - Brightness temperature of the forest for horizontal 
polarization and for a viewing angle of 50 degrees. 

Summary 

Examples of the brightness temperatures for the different classes 
and the effect of soil moisture and soil temperature are shown in 
Figu''es 17 thorugh 21. Note in Figures 18 and 19 that the brightness 

temperature of only three classes are dependent on soil moisture. 
These figures are for L-band and also demonstrate the effect of 
vegetation in decreasing the dependence on soil moisture. Figures 20 
and 21 show that at X-band there are only two classes that have a 
dependence on soil moisture. Also, the mixed bare and vegetative 

class is significantly less sensitive to sol moisture than the bare 
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FIGURE 21. Example of the dependence of brightness temperature on soil 
moisture for all classes at X-band vertical polarization. 



soil class. The effect of soil moisture on the brightness temperature 
of these classes at C-band lies between the effects shown at L-band 
and X-band. Figure 22 shows the effect of the soil temperature 
parameter on the brightness temperature of each class at X and L-band 
and \ 1 % soil moisture. 

SIMULATION RESULTS 
Simulation Description 

Rationale 

There are several scene and system parameters that can affect the 
ability to use a microwave radiometer for estimating soil moisture 
over extended scenes. This study was only concerned with the problems 
associated with the estimation of soil moisture assuming a spatially 
uniform soil moisture distribution over the scene under consider- 
ation. No consideration was given to the effects of soil moisture 
profile on the emission at the various microwave frequencies. This is 
a separate problem and is handled in another study. In addition, the 
spatial soil moisture distributions that naturally occur due to preci- 
pitation patterns and variations in soil properties were considered to 
impose separate restrictions on the resolution of an orbiting micro- 
wave radiometer and v/ere not considered in this study. However, these 
spatial variations are important and could very well be the limiting 
factor on resolution. Their effect can easily be considered using the 
simulation model. This should be investigated in a follow-on study. 

The scene and sensor parameters that were of concern to this 
study were scene heterogeneity and its relationshp to sensor resolu- 
tion, surface roughness, soil temperature, and sensor incident angle. 


62 


50 


0 10 

SOIL TEI 


FIGURE 22. Example of the dependent 
temperature parameter fc 





Specific simulation runs were made to provide data to address each of 
these factors. However, the most emphasis in this study centered 
around the effect of scene heterogeneity and resolution since these 
are the most critical parameters about which the least is currently 
known. 

The simulation of realistic scene geometry and composition were 
vital to the accuracy of the results of this study. As a result, the 
scene was simulated based on full frame Landsat images classified into 
six categories. These categories were water, bare soil, urban, mixed 
bare and vegetation, fully vegetated and forest. The description of 
how the original classifications were aggregated into these six, and 
the problems of classification consistency between Landsat frames is 
described in the section entitled Model Definition and Structure. 
These classes were chosen as being representati ve of the scene para- 
meters that are important in affecting the relationship between micro- 
wave emission and soil moisture. The pixel size over which the 
antenna integration was performed was 0,24 by 0.24 km. 

Test Runs 

In utilizing the simulated scene to siinulate microwave radiometer 
measurements for ana vses of the effects of scene heterogeneity, it 
was necessary to be careful in choosing the ground tracts of the 
radiometer flight path. The simulated scene contains very diverse 
ground cover from very dense forest in East Texas to sparse vegetation 
in Central Texas. Analysis of simulated radiometer measurements to 
determine the ability to estimate soil moisture from space can be 
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severely biased by the amount of vegetation contained in the radio- 
meter resolution elements. As a result, two ground tracts were chosen 
based on two genera] criteria. The first criteria was that the ground 
tract pass over areas of heavy forest vegetation as well as areas of 
sparser vegetation, and that the ground tracts pass over features that 
would be recognizable from the simulated radiometer measurement. One 
ground tract runs from just north of Waco, Texas southeastward to Lake 
Livingston, Texas. The other ground tract runs east and west from 
approximately Kerrville, Texas eastward to Houston, Texas and out into 
the Trinity Bay area. These ground tracts are shown on the urban 
class map contained in Appendix B. It will be shown below and can be 
seen in the class maps contained in Appendix B that these ground 
tracts cover areas that are predominantly vegetation. 

Numerous radioeter measurement simulations were computed for the 
two ground tracts described above. These simulations were run using 
the parameters documented in Table 4. The choice of parameters in 
Table 4 was based on the desire to determine the effects of resolution 
and frequency on the ability to estimate soil moisture with microwave 
radiometers over realistic and heterogeneous scenes. The simuation 
model computed both the vertical and horizontal brightness temperature 
for nadir angles from 0® to 50® in 10® increments. 

In order to quantify the effect of scene makeup on the microwave 
radiometer brightness temperature computation, the model was con- 
structed to compute and keep track of the percentage of each class 
contained in each radioemter footprint. These percentages are plotted 
as a function of range down the ground tract from west to east in 
Figures 23 through 28. Figures 23 and 24 contain the percentages of 
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TABLE 4. Simulation parameters used in Test Runs. 


Parameter 

Value 

Frequency 

L, C, X-band 

Soil Moisture 

0 %, 35 % 

Temperature Parameter 

10°C, 60°C 

Roughness Parameter 

0.3 

Antenna Footprint 

5 km, 20 !<m, 60 km 
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FIGURE 23. Percent of water, bare soil and forest classes in each antenna footprint as a function of 
range for the Waco to Livingston ground track at a resolution of 5 kilometers. 
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FIGURE 24. Percent of urban, mixed bare and vegetated and fully vegetated classes In each antenna 
footprint as a function of range for the Waco to Livingston ground track at a resolution 
of 5 kilometers. 
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FIGURE 25. Percent of water, bare soil and forest classes in each antenna footprint as a function of 
range for the Kerrville to Houston ground track at a resolution of 5 kilometers. 



ORIGINAL PAGE IS 
OF POOR QUALny 



.4 


70 


too. 130. 200. 2W. 300 . 330 . 400 . 430 . 300 . 830. 400. 

cni> 

FIGURE 26, Percent of urban, mixed bare and vegetated and fully vegetated classes In each antenna 
“ footprint as a function of range for the Kerrville to Houston ground track at a 5 

kilometer resolution. 
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for the Waco to Livingston ground track at a 20 kilometer resolution. 
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cent of water, bare soil and forest classes in each antenna footprint as a function of 
qe for the Waco to Livingston ground track at a 60 kilometer resolution. 




each of the classes in each antenna footprint for a 5 kilometer radio- 
meter resolution for the Waco to Livingston ground tract. Two water 
features show up distinctly on the plot in Figure 23. At a range of 
approximately 300 kilometers Lake Conroe shows up distinctly and at a 
range of approximately 525 kilometers Lake Livingston shows up 
distinctly. In addition, it can be seen that near Waco, the percent- 
age of forest in each footprint is on the order of 20% to 25%, while 
the percentage of forest begins to increase at approximately 375 kilo- 
meters of range and approaches 100% beyond 550 kilometers of range. 
Also in Figure 23, it can be seen that the amount of bare soil is very 
small, never going above approximately i0% of each footprint, except 
in the range of 325 kilometers to 374 kilometers where it averages 
approximately 25% of the footprint. In Figure 24, it can be seen that 
although the ground tract is over several cities, the percentage of 
urban area in each footprint is very small. It can also be seen in 
Figure 24 that the mixed bare and vegetated class and the fully vege- 
tated class made up a significant percentage of each footprint on the 
western end of the ground tract and decreases as a function of range. 

The same type of observations are made concerning Figures 25 and 
26 which show the percentage class of each 5 kilometer footprint for 
the Houston ground tract. Distinguishing features in this case are 
Trinity Bay which begins approximately 550 kilometers of range and the 
Houston urban area which spans the area from approximately 490 
kilometers to 550 kilometers of range. It can be seen that the 
Houston urban area does constitute approximately 40% to 50% of each 
footprint within the 500 to 550 kilometer range for the 5 kilometers 
antenna footprint. 
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Another factor of importance Is to note that as the antenna foot- 
print size is increased the percentage of each class making up each 
radiometer footprint remains approximately the same. Figures 27 and 
28 demonstrate this fact for an antenna footprint of 20 kilometers and 
60 kilometers for the Waco to Livingston ground tract and the water, 
bare soil and forest classes. It can be seen that the increased 
radiometer resolutions produce an averaging’ affect on the percent of 
class plots. It should also be pointed out that for a fixed ground 
tract length, the number of brightness temperature simulations 
decreases significantly as radiometer resolution is increased. 

Figure 23 through 28 illustrate the make-up of the simulated 
scene within the radiometer footprint spatially along the ground 
tracts. Figures 29 through 40 provide some additional statistical 
data concerning the make-up of the simulated scene along the ground 
tracts, but without the spatial information. These figures are fre- 
quency bar charts that show distribution of footprints along the 
ground tract in terms of the ground cover make up with the foot- 
prints. This distribution is computed in terms of 5% increments of 
total ground cover within each footprint. Figures 29 through 34 are 
for the Waco to Livingston ground tract, while Figures 35 through 40 
are for the Houston ground tract. All are computed for an antenna 
footprint of 5 kilometers. Again, it can be seen that by far the 
highest percentage of radioemter footprints is vegetation. The high- 
est percentage of bare soil in any footprint along the Houston ground 
tract is 25%, while the highest percentage of bare soil along the Waco 
to Livingston ground tract is 35%. In addition, by far the largest 
number of footprints contain 5% or less of bare soil. 
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FIGURE 29. Histogram of the percentage of the antenna footprint 

containing water for the Waco to Livingston ground track 
at a 5 kilometer resolution. 
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FIGURE 30. Histogram of the percentage of the antenna footprint 
containing bare soil for the Waco to Livingston ground 
track at a 5 kilometer resolution. 
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FIGURE 31. Histogram of the percentage of the antenna footprint 
containing mixed bare and vegetation for the Waco to 
Livingston ground track at a 5 kilometer resolution. 
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FIGURE 32. Histogram of the percentage of the antenna footprint 
containing full vegetation for the Waco to Livingston 
ground track at a 5 kilometer resolution. 
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FIGURE 33. Histogram of the percentage of the antenna footprint 

containing forest for the Waco to Livingston ground track 
at a 5 kilometer resolution. 
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urban for the Waco to Livingston ground track at a 5 kilometer 
resolution. 
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FIGURE 35. Histogram of the percentage of the antenna footprint containing 
water for the Waco to Livingston ground track at a 5 kilometer 
resolution. 
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Histogram of the percentage of the antenna footprint containing 
bare soil for the Kerrville to Houston ground track at a 5 
kilometer resolution. 
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Histogram of the percentage of the antenna footprint containing 
mixed bare and vegetation for the Kerrville to Houston ground 
track at a 5 kilometer resolution. 
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FIGURE 38. Histogram of the percentage of the antenna footprint containing 
full vegetation for the Kerrville to Houston ground track at a 
5 kilometer resolution. 


84 



Number of Antenna Footprints 


FKCQUENCy 

I 
I 
I 

45 ♦ 

I 
I 
I 
I 

40 * 

I 
I 
I 
i 

3S ^ 

I 
I 
I 
I 

30 * 

I 
I 
J 
I 

25 ♦ 

I 
I 
I 

20 ♦ 

I 
I 
I 
I 

IS ♦ 

I 
I 
I 
I 

10 ♦ *F 

I •« 

t *• 

I *• 

s ♦ 

I •• 

I ** 

I •« 

5 

FIGURE 39. 


0RJ61WAU PAGE IS 
OF POOR QUALITY 


(HOUSTON KUN 5KH RCSL) 
FRCGUENCV BAR CHART 


IB 

13134 HOtlOAV, JANUARY ZB, 1900 


* 


f 


10 



0 0 
00 

00 

00 

00 












0 0 
0 0 

0 0 
0 * 

* IF 

00 
0 0 
00 
0 0 


' 








• 


0 0 
0 0 

0 0 
00 

0 N 

Ik* 

00 

00 

• * 










0 0 

M 0 

• * 

0 0 

0 0 










0 0 

00 

** 

0 0 

0 0 










0 0 


• * 

00 

0 0 









m m 

0 0 

00 

* • 

00 


«* 








• • 


0 0 

* * 

00 

0 0 

*« 








** 

0 0 

0 0 

0 0 

00 

0 0 

• * 








m 

0 * 

00 

0 0 

0 0 

00 

** 








• 0 

0 0 

00 

0 0 

0 0 

0 0 

00 

** 







0 0 

00 

0 0 

00 

0 0 

00 


«« 







00 

00 


00 

0 0 

0 0 

00 

«« 







0 0 

0 0 



00 

0 0 

00 

** 







00 

0 0 

00 

00 

0 0 

0 0 

00 

«* 







00 

0 0 

0 k 

0 0 

0 0 

0 0 

00 

** 







0 0 

0 0 

0 0 

00 

0 0 

0 0 

00 

*• 







00 

00 


00 

00 

0 0 

00 


** 






00 

0 0 

0 0 

00 

00 

0 0 

00 

** 

** 






00 

0 0 

0 0 

00 


00 



• * 

00 

• • 



*• 

00 

00 

0 0 

00 

00 

00 

0 0 

0 0 

• * 

00 

• • 



00 

00 

0 0 

0 0 

00 

00 

0 0 

00 

00 


00 

*« 

• • 


00 

00 

00 

00 

00 

0 0 

0 0 

0 0 

00 

** 

0 0 

«« 

00 

00 

00 

15 

20 

ZS 

30 

35 

40 

45 

50 

55 

60 

65 

70 

75 80 85 

90 gs 100 


Percent of Footprint containing Forest 


Histogram of the percentage of the antenna footprint containing 
forest for the Kerrville to Houston ground track at a 5 
kilometer resolution. 
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FIGURE 40. Histogram of the percentage of the antenna footprint containing 
urban for the Kerrville to Houston ground track at a 5 
kilometer resolution. 
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Since the simulation model is reasonably expensive to run on the 
Texas ASM University Amdahl computers, only the two ground tracts were 
simulated in this study. However, in the interest of future analysis, 
all simulated outputs were recorded on magnetic tape and stored for 
future reference. 

Analysis 

Model Validity 

The steps in analyzing the simulated radiomter measurements were 
to verify that the simulation program behaves properly. It was impor- 
tant to demonstract that the brightness temperature computations were 
properly dependent on incident angle, temperature, surface roughness, 
soil moisture, and resolution. Figure 41 demonstrates the brightness 
temperature computation vs. incident angle at L-band for both horizon- 
tal and vertical polarizations from nadir to 50°. It can be seen in 
Figure 23 that the antenna footprint containing the largest percentage 
of bare soil occurred at a range of 350 kilometers in the Waco to 
Livingston ground track. In addition, it can also be seen in the same 
figure that the antenna footprint containing the largest percentage of 
water occurred at the ’range of 525 kilometers. Figure 41 demonstrates 
that the simulation model adequately predicts angular behavior of the 
brightness temperature in that the brightness temperature adequately 
responds to the scene makeup, within the antenna footprint as computed 
by the simulation model. 

Figure 42 demonstrates the effect of the temperature parameter 
computed by the simulation model. Figure 42 shows the horizontal 
brightness temperature computer for L-band at nadir for a 20 kilometer 
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INCIDENT ANGLE 


FIGURE 4i. Computation of brightness temperature as a function of incident 
angle for L-band and a 20 kilometer resolution for an antenna 
footprint containing predominantly bare soil and one containing 
predominantly water. 
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antenna footprint along the Waco to Livingston ground track for a tem- 
perature parameter of 10-C and a temperature parameter of 60®C. The 
effect of the 50® difference in temperature parameter is apparent. 
Also, the effect of two large water bodies are readily visible in the 
brightness temperature computation. Lake Conroe occurs at approxi- 
mately the 300 kilometer range and Lake Livingston occurs at approxi- 
mately the 525 kilometer range. It can also be seen that the effect 
of the temperature is different between footprints containing predomi- 
nantly forest and footprints with large percentages of bare soil. 
This is apparent by the difference in the two brightness temperature 
computations at a range of approximately 325 kilometers where the 
highest percent of bare soil occurs, and a range of approximate 550 
kilometers where the highest percentage of forest occurs. 

Figures 43 through 48 demonstrate the performance of the simula- 
tion model as a function of microwave frequency, antenna footprint 
size, and soil moisture. Figures 43, 44, and 45 are plots of horizon- 
tal brightness temperature computed at 35° incidence for an antenna 
footprint of 5 kilometers for L, C, and X-band respectively as a func- 
tion of range along the Waco to Livingston ground track for two values 
of soil moisture, 5% and 35%. These computations were made using a 
roughness factor of 0.3 which corresponds to a medium scale rough- 
ness. There are several observations which can be made concerning 
these three figures. First, in Figure 43, the effect of the increas- 
ing density of the forest from a range of approximately 350 kilometers 
to 600 kilometers is obvious. The large difference between the hori- 
zontal brightness temperature at 5% soil moisture and 35% soil mois- 
ture decreases as the forest cover density Increases. Beyond the 
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FIGURE 44, Horizontal brightness temperature computed at C-band and a 5 kilometer resolution for the 
Waco to Livingston ground track at two soil moistures, 5% and 35%. 
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FIGURE 45. Horizontal brightness temperature computed at X-band and a 5 kilometer resolution for the 
Waco to Livingston ground track at two soil moistures, 5% and 35%. 



FIGURE 46. Horizontal brightness temperature computed at L-band and a 20 kilometer resolution for 
the Waco to Livingston ground track at two soil moistures, 5% and 35%. 
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FIGURE 47. Horizontal brightness temperature computed at L-band and a 60 kilometer resolution for 
the Waco to Livingston ground track at two soil moistures, 5% and 35%. 
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range of 550 kilometers, where the percent of forest cover is in the 
90% range, the sensitivity to soil moisture is practically elimi- 
nated. By cotnparing Figure 43, 44, and 45 the effect of microwave 
frequency can easily be seen. It is obvious fron these figures that 
as the frequency goes higher, a difference between the brightness tem- 
perature at 5% soil moisture and 36% soil moisture decreases. This is 
due primarily to the effect of the vegetation cover that exists in the 
simulation scene. The largest percentages of bare soil occur around 
the range of 350 kilometers. Again, this Is obvious in Figures 43, 
44, and 45 since it is in this range interval that the largest differ- 
ence occurs between the 5% soil moisture computation and the 35% soil 
moisture computation for all microwave frequencies. 

The effect of antenna footprint size can be seen by comparing 
Figure 43, 46, and 47. The parameters used in generating these 
figures are identical except for the fact that the antenna footprint 
was increased from 5 kilometers to 20 kilometers and 60 kilometers, 
respecti vely . The effect of the largest footprint is obvious in the 
smoothing effect of the brightness temperature computations as a 
function of range. In addition, the ability to resolve physical 
features is diminished. In Figure 46, Lake Conroe and Lake Livingston 
are still visible, however, in Figure 47, Lake Conroe is practically 
unresol vable. However, the effect of the higher percentages of bare 
soil in the range around 350 kilometers still is visible in Figure 47 
for the 60 kilometer antenna footprint size. 

Figures 43 and 48 demonstrate the capability of the simulation 
model to compute horizontal and vertically polarized brightness 
temperatures. Figure 48 is identical to Figure 43 except that it pre- 
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sents verticatly polarized brightness temperature computations. It 
can be seen by comparing these two figures that the vertical polarized 
brightness temperature at both the 55 soil moisture and 355 soil mois- 
ture are several degrees higher than the horizontal brightness temper- 
ature. 

Analysis Approach 

As stated earlier, the objective of this study was to investigate 
the effects of scene heterogeneity on the ability of a microwave 
radiometer system to estimate soil moisture and to determine if the 
effects of scene heterogeneity are dependent upon the size of the 
antenna footprint. The most pleasing appraoch to accomplishing this 
objective is to quantify the effect of the percent ground cover of 
each class on the sensitivity of the brightness temperature computa- 
tion to soil moisture. Unfortunately , it is impossible to quantify 
the effects of each class independently of one another due to the di- 
mensionality of the problem. Therefore, the approach was to quantify 
the effects of each ground cover class on the sensitivity of the 
brightness temperature to soil moisture individually without regard to 
the other classes, and do this as a function of antenna footprint 
size. Care must be exercised in analyzing results presented in the 
manner. Obviously at low percentages of the ground cover class being 
investigated, there will be larger influences of other classes that 
have higher ground cover percentages. This will cause scatter in the 
results not due to the class being investigated. It is possible, 
however, to document the mean sensitivity to soil moisture for the two 
ground tracks considered and to approximate a quantitative description 
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of the effects of each Individual ground cover class on that sensitiv- 
ity. 

Analysis Results 

In this section the simulation results will be analyzed to demon- 
strate the effect of each scene class, microwave frequency, and 
antenna resolution on the sensitivity of brightness temperature to 
soil moisture. This will be done by first considering horizontal 
brightness temperature plotted as a function of percent class for two 
soil moistures, 5% and 35%. Figures 49 through 52 contain brightness 
temperature computation for the Waco to Livingston ground track at a 
35® incident angle and a 5 kilometer resolution. These plots demon- 
strate the scatter In the brightness temperature computations when the 
computations for each ground resolution cell are viewed Independently. 

Figure 49 shows L-band horizontal brightness temperatures for 
bare soil plotted as a function of percent class of bare soil. The 
general effect of the soil moisture difference is obvious. However, 
between 0% and 8% bare soil there are some very low brightness temper- 
atures. These particular footprints are those that occurred over Lake 
Conroe and Lake Livingston, thus these low brightness temperatures are 
due to water. Note also, that the largest amounts of scatter In the 
brightness temperatures occur In the low percentage range of bare 
soil, and that as the percentage of bare soil Increases the scatter In 
the brightness temperature computation decreases. In addition, the 
largest scatter occurs for the brightness temperature computations at 
35% soil moisture. These effects are due to the fact that at low 
percentages of bare soil the primary class consitutents within the 
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FIGURE 49. Horizontal brightness temperature computed at L-band for the Waco to Livingston ground 
track at a 5 kilometer resolution plotted as a function of percent class of bare soil at 
two soil moistures, 5% and 35%. 
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FIGURE 50. Horizontal brightness temperature computed at L-band for the Waco to Livingston ground 
track at a 5 kilometer resolution plotted as a function of percent class of mixed bare 
and vegetation at two soil moistures, 5S and! 35S. 
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FIGURE 51. Horizontal brightness temperature computed at L-band for the Waco to Livingston ground 
track at a 5 kilometer resolution plotted as a function of percent class of fully 
vegetated at two soil moistures, 5% and 35%. 
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resolution footprint is some other class. No information is pro- 
vided on which class in this particular plot. However, by looking at 
Figures 23 and 24 that show the percentage class distribution as a 
function of range down the Waco to Livingston ground track, it is 
obvious that the majority of the footprint for the low percentage of 
bare soil resolution are some fonn of vegetation, primarily, fully 
vegetated and forested. 

Figure 50 shows the same type of information except for the mixed 
vegetation class. It can be seen that above approximately 20% mixed 
bare soil and vegetation the scatter in the brightness temperature 
computations are significantly decreased. Figure 51 shows the 
brightness temperatures for 5% and 35% soil moisture plotted as a 
function of the fully vegetated class. Note that there are antenna 
footprints where the percentage of this class is much higher than that 
for bare soil and the mixed vegetation class. Also, the scatter for 
the 35% soil moisture brightness temperature computation continues to 
much higher percentages of the fully vegetated class than previously 
seen for the mixed vegetation class or the bare soil class. This is 
due to the fact that the resolution elements that contained less than 
approximately 40% of the fully vegetated class contain very high 
percentages of the forest class, usually above 50% as can be seen in 
Figures 23 and 24. 

Figure 52 shows the same type of plot for the forest class. The 
effect of increasing the percentage forest in each resolution element 
is obvious. It can be seen that above approximately 35% to 40% of the 
forest class there is severe degradation in the ability to distinguish 
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between the two soil moisture contents. In addition, below approxi- 
mately the 40% point there are some extremely low brightness tempera- 
ture computations. Again, this Is due primarily to Lake Livingston 
and Lake Conroe. 

Plots such as those shown in Figure 49 through 52 will not be 
provided for the urban and water class since there was not a good 
enough distribution of resolution elements with different percent 
classes to present the data In this manner. 

Figures 49 through 52 show brightness temperature computations 
for L-band at 5 kilometers for each of the four major classes. 
Figures 53 and 54 show the same type computations for the class mixed 
vegetation at C-band and X-band, respectively. In Figure 53 It can be 
seen that the same general comments that were made concerning Figure 
50 (which showed mixed vegetation class for L-band) can be made. The 
only difference is that the magnitude difference between the horizon- 
tal brightness temperature computed at 5% soil moisture and the 
brightness temperature computed at 35% soil moisture is less at C-band 
than it was at L-band. In addition, the percent of the mixed vegeta- 
tion class must get above approximately 20% In order for the scatter 
in the brightness temperature for the 35% soil moisture to decrease 
significantly. In Figure 54, which is the same plot for X-band, it 
can be seen that the effect of the vegetation is more severe. In 
fact, the percent class of the mixed vegetation must be above approxi- 
mately 35% for the scatter to significantly decrease* It should also 
be noted that the ability to discriminate soil moisture increases as 
the percent of the mixed vegetation class increases at all three 
frequencies. This is due to the fact that the forest class predomi- 
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FIGURE 54. Horizontal brightness temperature confuted at X«band for the Waco to Livingston ground 
track at a 5 kilometer resolution plotted as a function of percent mixed bare and 
vegetation class at two soil moistures, S% and 35%. 



nates the jnake of the resolution elements for the low percentages of 
the mixed vegetation class. 

Figures 55 and 56 will be used to demonstrate the effect of 
resolution on the scatter seen In the brightness temperature computa- 
tions. Figure 55 shows the horizontal brightness temperature computed 
for L-band at 5X and 355t soil moisture for the mixed vegetation 
class and a 20 kilometer resolution. This figure can be compared to 
Figure 50 which Is the same plot for the 5 kilometer resolution size. 
It can be seen that the general trends In both figures are the same 
although the scatter in the data for the 20 kilometer footprint Is 
significantly less. Although there Is less scatter, there Is a cor- 
responding fewer number of points In Figure 55 than In Figure 50. 
This demonstrates the averaging effect of simply increasing the 
resolution size of the antenna. Increasing resolution is very nearly 
equivalent to averaging measurements taken at smaller resolution. The 
effect of going to a 60 kilometer resolution element Is essentially 
the same, additional averaging and fewer points. This result indi- 
cates that in general the ground cover classes behave as If they are 
randomly oriented over the scene. Figure 56 demonstrates the s^mie 
phenomena as Figure 55. The four points at the 45?i forest class cor- 
respond to calculations made over Lake Livingston as can be seen from 
Figure 27. 

Figures 49 through 56 provide a qualitative indication of the 
distribution of the brightness temperature computations as a function 
of scene class. It is Informative to investigate the distribution of 
the brightness temperature computations for 5/t soil moisture and 35X 
soil moisture for the entire scene without regard to Individual 
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FIGURE 55. Horizontal brightness ter^ierature computed at L-band for the Waco to Livingston ground 
track at a 20 kilometer resolution plotted as a function of percent mixed bare and 
vegetation class at two soil moistures, 5% and 35%. 
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FIGURE 56. Horizontal brightness temperature computed at L-band for the Waco to Livingston ground 
track at a 20 kilometer resolution plotted as a function of percent forest class at two 
soil moistures, 5% and 35%. 




classes. Figures 57 and 58 are frequency histograms for the bright- 
ness temperature computed at L-band at an Incident angle of 35° for 
the Waco to Livingston ground track at horizontal polarization. 
Figure 57 was computed for a 5 kilometer ground resolution. It can be 
seen that the distribution of brightness temperatures for a S% soil 
moisture tightly clusters between 250°K and 275°K. The brightness 
temperatures for the 35X soil moisture has a peak at 200°K but spreads 
up to 275°K. It will be seen later that this spreading effect causes 
the confidence Interval on the sensitivity of the brightness tempera- 
ture to soil moisture to be larger than that for the 60 kilometer 
resolution. Figure 58 shows the same frequency histogram for the 60 
kilometer resolution. Again, for the 5.i soil moisture the brightness 

temperature cluster tightly between 250 and 275°K. However, for the 

/ 

355i soil moisture, brightness temperature is a multimodal 
distribution. One peak occurs at 250°K while the other occurs at 
approximately 230°K. 

Up to now we have only qualitatively investigated the dependence 
of brightness temperature on soil moisture by analyzing brightness 
temperature computations for 5% soil moisture and 35% soil moisture 
over the same ground track. One method of quantifying the sensitivity 
of the brightness temperature to soil moisture is to compute the slope 
of the best fit straight line between the brightness temperatures com- 
puted for 5% soil moisture and 35% soil moisture. Since It is known 
from other investigations that the microwave brightness temperature Is 
linearly related to soil moisture, this slope can be termed the 
sensitivity of the brightness temperature to soil moisture in ®K per 
percent soil moisture. Figures 59 through 61 are plots of the sensl- 
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tivity of the brightness temperature to soil moisture as a function of 
percent class of bare soil for the 20 kilometer resolution along the 
Waco to Livingston ground track computed at horiiontal polarization 
for L-band, C-band. and X-band respectively. It can seen that there 
is considerable scatter in the sensitivity up to approximately 20^ 
bare soil. Above approximately 20X bare Soil the sensitivity at 
L-band is approximately 1.75®K for percent soil moisture. Again it 
can be noted that the scatter in the sensitivity below 20% bare soil 
is due to a very high percentage of forest and fully vegetated in the 
ground resolution elements (individual antenna footprints). This can 
be verified from Figures 23 and 24. In Figure 23 it can be seen that 
only between the range of 325 kilometers and 375 kilometers does the 
bare soil go above approximately 20% of the antenna footprint. Also 
in Figure 23 it can be seen that above a range of approximately 375 
kilometers the majority of the footprint is forest and in Figure 24 
below a range of approximatey 325 kilometer the majority of the foot- 
print is fully vegetated. The same plot for C-band is shown in Figure 
60. The same effects are seen to exist, however, the sensitivity of 
the brightness temperature to soil moisture above approximately 20% 
bare soil only reaches approximately l,25°K/percent soil moisture. 
This is due to the stronger effect of the vegetation at C-band than at 
L-band. Figure 61 shows the sante information for X-band. Again, the 
same effects are seen to exist, but effects of other classes besides 
bare soil are obviously more severe. There is more scatter in the 
sensitivity computations than for L-band and C-band above 20% bare 
soil. The average sensitivity at 30% bare soil is approximately 
l.l*’K/percent soil moisture. 
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Figures 59 through 61 demonstrated the effect of microwave fre- 
quency on the sensitivities to soil moisture using the bare soil class 
as an example. Figure 62 through 64 will be used to demonstrate the 
effect of the classes on the sensitivity to bare soil using L^-band 
horizontal polarization, 20 kilometer footprint size and the Waco to 
Livingston ground track. Figure 62 shows the sensitivity plotted as 
the function of percent class for L-band and the mixed vegetation 
class. By comparing this figure to Figure 59 it can be seen that the 
results are nearly identical to bare soil except that the percentage 
of the mixed vegetation class within the footprint goes up to nearly 
40^6 while the percentage of bare soil only went slightly above 30%. 
In addition, the sensitivity climbs to approximately 2°K/percent soil 
moisture. This sensitivity is fairly large since the 40% mixed 

vegetation class occurs between the ranges of 125 kilometers to 140 
kilometers and a 150 kilometers to 175 kilometers as seen in Figures 

23 and 24. In these range intervals the percentage of forest is 

somewhere between 20% and 25%. 

Figure 63 shows the same information except for the fully 

vegetated class. The character of this plot looks slightly different 
than the previous ones in that it appears that the sensitivity 
linearly climbs as a function of percent class although there is some 
scatter between 38% and 48% class. This is due to the fact that below 
38% class of fully vegetated soil the primary other component of the 
ground resolution footprint is forest as can be seen from Figure 23 
and 24. This occurs at a range of above 350 kilometers where the 
percent of forest ranges anywhere from 50 to 90%. The scatter in 
Figure 63 toward higher sensitivities between 38% class and 48% class 
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FIGURE 63. Soil moisture sensitivity of the L-band horizontally polarized brightness temperature for 
the Waco to Livingston ground track at a 20 kilometer resolution plotted as a function of 
percent fully vegetated class. 
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occurs due to the significant percentage of bare soil within those 
resolution footprints. The 38X to 48% of fully vegetated soil occurs 
between the ranges of approximately 300 to 350 kilometers where the 
bare soil approaches its maximum. At any rate* it can be seen that 
since the sensitivity above 50% of the fully vegetated class is 
greater than 1.5®K percent soil moisture, the fully vegetated class is 
not the controlling effect on the sensitivity of brightness 
temperature to soil moisture at least for L-band. 

Figure 64 shows the same information for the forest class. In 
the models used to compute the brightness temperature for each class, 
it was assumed that the forest class had no sensitivity to soil mois- 
ture. As a result Figure 64 is an expected result. As the percent of 
forest increases the sensitivity to soil moisture linearly decreases 
from approximately 2°K/percent soil moisture at 15% of the forest 
class to zero at 100% of the forest class. There were so few 
resolution elements with a significant percentage of urban or water 
class that it was pointless to attempt to make similar plots for those 
classes. The general conclusion that can be drawn from Figures 59 
through 64 is that the forest class is probably the controlling effect 
on the sensitivity of the brightness temperture to soil moisture. 

Figures 59 through 64 were used to show the sensitivity of the 
brightness temperature to soil moisture as a function of the percent- 
ages of each scene constituent within an antenna footprint. It is 
also instructive to compare the sensitivity of the brightness 
temperature to soil moisture for all ground resolution elements 
Without regard to the constituency of the scene within each antenna 
footprint as a function of frequency and resolution. Figure 65 is a 
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plot of the sensitivity of the horizontally polarized brightness 
temperature to soil moisture as a function of resolution and 
frequency without regard to the class make up of the scene. 

Since It was previously shown that the forest class is the con- 
trolling effect on the sensitivity of the brightness temperature to 
soil moisture, there Is obviously significant effects of the high per- 
centages of forest with the simulated scene on the average sensitivity 
to soil moisture presented In Figure 65. In an operational satellite 
system put up for the purpose of estimating soil moisture using a 
microwave radiometer, the location of each antenna footprint should be 
known. It would not be difficult to map the areas of significant 
forest cover such that footprints containing percentages of forest 
cover over a certain threshold are neglected. 

The improvement in sensitivity to soil moisture when it is pos- 
sible to partition ground resolution elements by the percentage of 
forest within each element is investigated in Figures 66 and 67. In 
order to obtain enough points to be statistically significant, both 
the Waco to Livingston and Houston ground tracks were used to compute 
the numbers shown in Figures 66 and 67. Figure 66 shows a plot of 
sensitivity in per percent soil moisture plotted as a function of 
resolution and frequency. The sensitivity plotted in Figure 66 is an 
average computed only for antenna footprints containing less than 40% 
forest cover. In addition, 95% confidence intervals for these average 
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sensitivities are also provided. It can be seen that at Uband the 
sensitivity is approximately 1.5®K/percent soil moisture for all 
resolution with a very slight increase as the resolution increases. 
At C-band the sensitivity is approximately 0.8°K/percent soil moisture 
and at X-band approximately 0.5®K/percent soil moisture. Figure 68 
provides the same information except at the sensitivities plotted in 
figure 68 were computed only for antenna footprints containing less 
than 30 percent forest cover. A major effect is that the sensitivity 
at L-band increases from approximately 1.5®K/percent soil moisture to 
1.64®K/percent soil moisture fran the 5 kilometer resolution to the 60 
kilometer reslution. Similar effects occur at C-band and X-band but 
to a lesser degree. In both Figure 66 and 67 the tremendous sensitiv- 
ity reduction at C-band and X-band relative to L-band is due to the 
effects of vegetation. Another very important fact to notice is that 
the confidence interval computed for the average sensitivity is the 
largest at the 20 kilometer resolution. For all frequencies the con- 
fidence intervals are comparable at the 5 kilometer and 60 kilometer 
resolutions. There were not enough data points for antenna footprints 
containing only lOJi and 20% forest cover to provide similar plots 
based on these partitioning percentages. 

CONCLUSIONS 

A simulation model of an orbiting microwave radiometer has been 
implemented and demonstrated. A realistic scene was also constructed 
over which the radiometer could be arbitrarily flown in a line scan or 
side to side scanning mode* The ground scene is based on classified 
Landsat images and thereby provides realistic ground classes os'" well 
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as geometries. It was shown that the realism of the scene make up In 
both ground cover class and geometry is critical to the accuracy of 
the simulation results. The model is capable of computing brightness 
temperatures for L-band, C-band and X-band frequencies as a function 
of orbit and antenna characteristics. 

An analysis of the effects of scene heterogeneity and antenna 
resolution on the sensitivity to soil moisture was performed. It was 
assumed in executing this analysis that the soil moisture was uniform 
over the entire scene (i.e. no precipitation patterns were overlaid on 
the scene) such that effects of the different scene components on the 
sensitivity to soil moisture could be Investigated. A sensitivity to 
soil moisture as a function of the percentage of each ground cover 
class was shown. It was demonstrated that the forest cover class was 
the limiting factor on the sensitivity of brightness temperature to 
soil moisture. The average sensitivity to soil moisture as a function 
of frequency and resolution without regard to ground cover classes was 
also shown. In addition, since the forest cover class was the limit- 
ing factor on sensitivity to soil moisture, average sensitivities for 
each frequency and resolution were computed based on partitioning out 
antenna footprints containing greater than 40% of the forest cover 
class. It was shown under this condition that the L-band frequency 
could achieve a resolution of 1.5®K/percent soil moisture, C-band a 
sensitivity of 0.8®K/percent soil moisture and X-band a sensitivity of 
0.5®K/percent soil moisture. The significant reduction as frequency 
increases is due to the effects of vegetation. It should be pointed 
out that the scene was made up primarily of the fully vegetated and 
forest cover classes. The largest percentage of the bare soil class 
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was in the range of 305i to 35% for only approximately 10% of the 
ground track simulated. 

Another important result was that the effect of increasing reso- 
lution size (antenna footprint size) is to provide a data averaging 
effect. However, it was seen that increasing resolution did improve 
the average sensitivity. This is most likely due to the manner in 
which the larger antenna footprint size averages the forest and fully 
vegetated classes. In addition, it was seen that the confidence 
interval on the average sensitivity computation was the largest for 
the 20 kilometer footprint and smallest for the 5 kilometer and 60 
kilometer footprints, which were approximately the same. 

This model is a very useful tool for a very wide ranging set of 
investigations. Although atmospheric effects were not considered in 
this study, the model is well suited for including an atmospheric 
model. In addition, the model would be well suited for quantitatively 
testing soil moisture estimation algorithms derived from other 
studies. It cart also be used to develop test efficient soil moisture 
estimators based on the Kalmen filtering approach. 
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APPENDIX A 


TEXAS PARKS & WILDLIFE MAPTAP VALUES 
VS 

TAMU CLASS ASSIGNMENTS 


A limited number of vegetation type maps were produced in color. 
Additional information contained on the maps may be obtained by 
contacting the authors of this report. 
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torsicana i/w.i6ZZ5 


legend Nomenclature and Color Assignments 
for Respective Values on MAPTAP 0062 


MAPTAP 

Value 

Legend Nomenclature 

TAMU 

CLASS 

0 

Unclassified 


1 

Grasses 

5 

2 

Grasses 

5 

3 

Grasses 

5 

4 

Post Oak-Black Hickory Forest 

6 

5 

Mesquite-Elm Parks 

6 

6 

Post Oak-Black Hickory Forest 

6 

7 

Mesquite Woods 

6 

8 

Pecan-Elm/Water Oak-Elm/Elm-'Hackberry Forest 

6 

9 

Pecan-Elm/Water Oak-Elm/Elm-Hackberry Forest 

6 

10 

Pecan-Elm/Water Oak-Elm/Elm-Hackberry Forest 

6 

11 

Pecan-Elm/Water Oak-Elm/Elm-Hackberry Forest 

6 

12 

Loblolly Pine-Sweetgum Forest 

6 

13 

Post Oak-Black Hickory Forest 

6 

14 

Post Oak-Black Hickory Forest 

6 

15 

Post Oak-Black Hickory Forest 

6 

16 

Loblolly Pine-Sweetgum Forest 

6 

17 

Water 

1 

18 

Water 

1 

19 

Water 

1 

20 

Water 

1 

21 

Crops 

4 
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Corsicana (2) 


fiAPTAP TAMU 

Value Legend Nomenclature CLASS 


22 

Crops 

4 

23 

Crops 

4 

24 

Crops 

4 

25 

Crops 

4 

26 

Crops 

4 

27 

Crops 

4 

28 

Crops 

4 

29 

Sparsely Vegetated/Urban 

2 

30 

Sparsely Vegetated/Urban 

2 

31 

Sparsely Vegetated/Urban 

2 

32 

Sparsely Vegetated/Urban 

2 

33 

Sparsely Vegetated/Urban 

2 

34 

Sparsely Vegetated/Urban 

2 

35 

Sparsely Vegetated/Urban 

2 

36 

Sparsely Vegetated/Urban 

2 

126 

Latitude-Longitude Tick Marks 


127 

Background 
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Kerrvine 1706-16343 


Legend Nomenclature and Color Assignments 
for Respective Values on MAPTAP 0068 

MAPTAP * TAMU 

Value . Legend Nomenclature CLASS 

0 Unclassified 

1 Grasses 5 

2 Grasses 5 

3 Grasses 5 

4 Live Oak-Ashe Juniper/Llve Oak-Mesquite Parks 5 

(Sparse) 

5 Sparsely Vegetated/Urban 2 

6 Live Oak-Mesquite Park (Sparse) 5 

7 Live Oak-Ashe Juniper/Live Oak-Mesquite Parks 5 

(Sparse) 

8 Live Oak-Ashe Juniper/Live Osk-Texas Oak/Live 6 

Oak-Mesquite Parks (Dense) 

9 Live Oak-Ashe Juniper/Live Oak-Texas Oak/Live 6 

Oak-Mesquite Parks (Dense) 

10 Live Oak-Ashe Juniper/Live Oak-Texas Oak/Live 6 

Oak-Mesquite Parks (Dense) 

11 Live Oak-Ashe Juniper Woods 6 

12 Live Oak-Ashe Juniper Woods 6 

13 Live Oak-Ashe Juniper Woods 6 

14 Pecan-Elm Forest 6 

15 Pecan-Elm Forest 6 

16 Live Oak-Ashe Juniper/Live Oak-Texas Oak/Live 6 

Oak-Mesquite Parks (Dense) 

17 Crops 4 

18 •’ Crops 4 

19 Crops 4 
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MAfTAP TAMU 

Value legend Nomenclature CLASS 

2G Water 1 

21 Sparsely Vegetated/Urban 2 

22 Sparsely Vegetated/Urban 2 

23 Cloud Cover 

24 Cloud Cover 


To be designated by specified colors but not included in the legend 

126 Latitude-Longitude Tick Marks 

127 Background 



Brownvraod 1706-16341 


Class Arrangement as Desired on Map Legend 
(MAPTAP No. 0059) 

NOTE: Class names and color assignment now occur only once on the 
MAPTAP and follow the Legend. 


MAPTAP TAMU 

Value Legend Nomenclature CLASS 

1 Grasses 5 

2 Oak/Mesquite/ Juniper/Mixed Parks (Sparse) 5 

3 Oak/Mesquite/Juniper/Mixed Parks (Dense) 6 

4 Oak/Mesquite/iluniper/Mixed Woods 6 

5 Crops 4 

6 Water 1 

7 Sparsely Vegetated/Urban/Crops 2,3 

8 Unclassified 


To be designated by colors but not included in the Legend. 

126 (Latitude-Longitude Tick Marks) 

127 (Background) 


Vegetation Type Map 







Austin 1363-16362 


Legend Nomenclature and Color Assignments 
fo>^ Respective Values on MAPTAP 0060 


MAPTAP 

Value 

Legend Nomenclature 

TAMU 

CLASS 

0 

Unclassified 


1 

Grasses 

5 

2 

Grasses 

5 

3 

Grasses 

5 

4 

Live Oak-Ashe Juniper/Post Oak-Live Oak Parks 
(Dense) 

6 

5 

Live Oak-Ashe Juniper/Post Oak-Live Oak Elm- 
Hackberry Parks (Sparse)/Grasses 

5 

6 

Live Oak-Ashe Juniper/Post Oak-Live Oak/El m- 
Hackberry Parks (Sparse)/Grasses 

5 

7 

Live Oak-Mesquite Parks (Sparse) 

5 

8 

Live Oak-Ashe Juniper/Post Oak-Live Oak/Elm- 
Hackberry Parks (Sparse)/6rasses 

5 

9 

Ashe Juniper Parks (Dense) 

6 

10 

Live Oak-Ashe Juniper/Post Oak-Live Oak/Elm- 
Hackberry Parks (Sparse )/6rasses 

4 

11 

Live Oak-Ashe Juniper/Post Oak-Live Oak/Eln- 
Hackberry Parks (Sparse)/Grasses 

4 

12 

Live Oak-Ashe Juniper Woods 

6 

13 

Post Oak-Blackjack Oak/El m-Hackberry Woods 

6 

14 

Pecan-Elm/Live Oak Forest 

6 

15 

Crops 

4 

16 

Crops 

4 

17 

Crops 

4 

18 

Crops 

4 

19 

Sparsely Vegetated/Urban 

nr 

2 
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MAPTAP TAMIJ 


Value 

Legend Nomenclature 

CLASS 

20 

Crops 

4 

21 

Crops 

4 

22 

Live Oak-Ashe Juniper Woods 

6 

23 

Cloud Cover 

0 

24 

Cloud Cover 

0 

25 

Water 

1 

26 

Water 

1 

27 

Water 

1 

28 

Sparsely Vegetated/Urban 

2 

29 

Sparsely Vegetated/Urban 

2 

30 

Sparsely Vegetated/Urban 

2 

31 

Sparsely Vegetated/Urban 

2 

32 

Sparsely Vegetated/Urban 

2 

33 

Post Oak-Eastern Redcedar Woods 

6 

34 

Post Oak-Eastern Redcedar Parks (Dense) 

6 

To be designated by specified colors but not included 

in legend. 

126 

Latitude-Longitude Tick Marks 


127 

Background 
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Waco 1759-16263 


Class Arrangement as Desired on 
Map Legend (MAPTAP No. 0086) 


NOTE: Class names and color assignnment now occur only once on the 
MAPTAP and follow the Legend. 


MAPTAP TAMU 

Value Legend Nomenclature CLASS 

1 Grasses 5 

2 Oak/Mesquite/ Juniper/Mixed Parks (Sparse)/ 5 

Grasses 

3 Oak/Mesquite/Juniper/Mixed Parks (Dense) 6 

4 Oak/Mesquite/Juniper/Mixed Woods 6 

5 Pecan-Elm/Live Oak Forest 6 

6 Crops 4 

7 Water 1 

8 Sparsely Vegetated/Urban 2,3 

9 Cloud Cover 

10 Unclassified 


To be designated by colors but not included in the Legend. 

126 (Latitude-Longitude Tick Marks) 

127 (Background) 






Lufkin 1757-16151 


Legend Nomenclature and Color Assignments 
for Respective Values on mAPTAP 5o4^^ 

MAPTAP TAMU 

Value Legend Nomenclature CLASS 

0 Unclassified 0 

1 Grasses 5 

2 Grasses 5 

3 Water Oak-Elm Forest 6 

4 Shortleaf Pine-Post Oak/Loblolly Pine- 6 

Sweetgum/Post Oak-Black Hickory Forest 

5 Shortleaf Pine-Post Oak/Loblolly Pine- 

Sweetgum/Post Oak-Black Hickory Forest 6 

6 Shortleaf Pine-Post Oak/Loblolly Pine- 6 

Sweetgum/Post Oak-Black Hickory Forest 

7 Loblolly Pine-Shortleaf Pine Forest 6 

8 Loblolly Pine-Shortleaf Pine Forest 6 

9 Loblolly Pine-Slash Pine Young Forest/ 

Eastern Mixed Hardwood Brush/Woods 6 

10 Water 1 

11 Sparsely Vegetated/Urban 2,3 

12 Crops 4 

13 Crops 4 

14 Crops 4 

To be designated by specified colors but not included in the legend. 

126 Latitude-Longitude Tick Marks 

127 Background 


-<F 
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Houston 1703-16173 


Legend Nomenclature and Color Assignments 
for Respective Values on MAPTAP 0031 


Legend Nomenclature 


Unclassified 0 

Grasses 5 

Grasses 5 

Grasses 5 

V/ater Oak-Elm/Pecan-Elm/Willow Oak-Blackgum 6 

Forest 

Post Oak-Black Hickory Forest 6 

Post Oak-Black Hickory Forest 6 

Loblolly Pine-Sweetgum/Loblolly Pine- 6 

Shortleaf Pine Forest 

Loblolly Pine-Sweetgum/Loblolly Pine- 6 

Shortleaf Pine Forest 

Loblolly Pine/Sweetgum/Loblolly Pine- 6 

Shortleaf Pine Forest 

Water Oak-Elm/Pecan-Elm/Willow Oak- 6 

Blackgum Forest 

Loblolly Pine-Sweetgum/Loblolly Pine- 6 

Shortleaf Pine Forest 

Water Oak-Elm/Pecan-Elm/Willow Oak- 6 

Blackgum Forest 

Loblolly Pine-Slash Pine Young Forest/ 6 

Eastern Mixed Hardwood Brush/Woods 

Loblolly Pine-Sweetgum/Loblolly Pine- 6 

Shortleaf Pine Forest 

Crops 4 

Crops 4 

Crops 4 


Houston (2) 


MAPTAP TAMU 


Value 

Leqend Nomenclature 

CLASS 

18 

Crops 

4 

19 

Crops 

4 

20 

Crops 

4 

21 

Crops 

4 

22 

Cloud Cover 

0 

23 

Elm-Hackberry Woods/Baccharis Brush 

5 

24 

Elm-Hackberry Woods/Baccharis Brush 

5 

25 

Urban/Sparsely Vegetated 

2 

26 

Urban/Sparsely Vegetated 

2 

27 

Urban/Sparsely Vegetated 

2 

28 

Urban/Sparsely Vegetated 

2 

29 

Urban/Sparsely Vegetated 

2 

30 

Urban/Sparsely Vegetated 

2 

31 

Urban/Sparsely Vegetated 

2 

32 

Urban/Sparsely Vegetated 

2 

33 

Urban/Sparsely Vegetated 

2 

34 

Marsh/Cultivated Wetlands 

1 

35 

Marsh/Cultivated Wetlands 

1 

36 

Loblolly Pi ne-Sweetgum/ Loblolly Pine- 
Shortleaf Pine Forest 

6 

37 

Water 

1 

38 

Loblolly Pine-Sweetgum/Loblolly Pine- 
Shortleaf Pine Forest 

6 

39 

Urban/Sparsely Vegetated 

2 

40 

Cloud Cover 

0 

41 

Cloud Cover 

0 
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Houston (3) 


To be designated by specified colors but not included in the legend. 

126 Latitude-Longitude Tick Marks 

127 Background 
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Bryan 1308-16311 


Legend Nomenclature and Color Assignments 
for Respective Values on MAPTAP 0067 


MAPTAP 

Value 

Legend Nomenclature 

TAMU 

CLASS 

0 

Unclassified 

0 

1 

6rasses/Elm-Hackberry Parks 

5 

2 

Post Oak-Blackjack Oak/Elm-Hackberry 
Woods/Forest 

6 

3 

Grasses/Elm-Hackberry Parks 

5 

4 

Grasses/Elm-Hackberry Parks 

5 

5 

Grasses/Elm-Hackberry Parks 

5 

6 

Grasses/Elm-Hackberry Parks 

5 

7 

Post Oak-Blackjack Oak/Elm-Hackberry 
Woods/Forest 

6 

8 

Sparsely Vegetated/Urban 

2 

9 

Crops 

4 

10 

Post Oak-Blackjack Oak/Elm-Hackberry 
Woods/ Forest 

6 

11 

Water Oak-Elm/Pecan-Elm Forest 

6 

12 

Post Oak-Blackjack Oak/Elm-Hackberry 
Woods/Forest 

6 

13 

Crops 

4 

14 

Grasses/Elm-Hackberry Parks 

5 

15 

Crops 

4 

16 

Water Oak-Elm/Pecan-Elm Forest 

6 

17 

Loblolly Pine-Pi nst Oak Forest 

6 

18 

Crops 

4 

19 

Crops 

4 
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Bryan (2) 


MAPTAP 

Value 

Lenend Nomenclature 

TAMU 

CLASS 

20 

Water Oak-Elm/Pecan-Elm Forest 

6 

21 

Sparsely Vegetated/Urban 

2 

22 

Sparsely Vegetated/Urban 

2 

23 

Sparsely Vegetated/Urban 

2 

24 

Water 

1 

25 

Loblolly Pine-Post Oak Forest 

6 

26 

Water 

1 

27 

Crops 

4 

28 

Loblolly Pine-Post Oak Forest 

6 

29 

Crops 

4 

30 

Crops 

4 

31 

Crops 

4 

32 

Grasses/Elm-Hackberry Parks 

5 

33 

Sparsely Vegetated/Urban 

2 

34 

Sparsely Vegetated/Urban 

2 

35 

Crops 

4 

36 

Crops 

4 

37 

Crops 

4 

38 

Sparsely Vegetated/Urban 

2 

39 

Crops 

4 


To be designated by specified colors but not included in the legend. 

126 Latitude-Longitude Tick Marks 

127 Background 
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APPENDIX B 


FULL SCENE GREY-SCALE riAPS BY CLASS 


original page is 

OF POOR QUALITY 



Figure B-1. Image of ground scene high-lighting pixels classified as "open water 
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Figure B-6. Grey-scale map of pixels classified as "fully vegetated 








APPENDIX C 


COMPUTATION FLOW OF COMPUTER ALGORITHM 


^ 7tart ^ 



fiNITIALIZED 
VARIABLES 


<!>, 0. B» 

®no* '^0* 

Xo =825.. Yo=0.. Zo-0 
AX = AY = 1., AZ = 0. 
Aa = Ay = 0., A3 = 1. 
A6 =l./57.3, C = 273. 
N = 2496, M = 1650, 

S = 0.24 


Rso“ Htdn 9 no 
“0 " Xo-RsoCosyo 

3o = Yo"Rso^^ •^'^0 
YO = H 


i 


OPEN 
UNIT 8 


0 


EVALUATE INITIAL NADIR COORDINATES 


KEY 


$, azimuth, - 
0 , roll , -90' 


■90° _< <1* £ 90° 

’ £ 0£ 90° 

B, beam width, 1° £ B £ 30° 
9f,o» incidence, 0° £0no 1 50° 

H'o* radial, 0° £ < 359° 

H, altitude, 50 £yq £500 KM 
N, scene length, 2496 units 
M, scene width, 1650 units 

C, background temp., 273°K 

S, scale factor, 0.24 KM/unit 
II, nadir step, l.£ll ^16 
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RESET NU,DE, 
OMX*0MN’#1X* 

4>mn» Rs- 

STEP a,g,Y 



SET LIMITS 
OF INTEG. 
II. 12 
K1.K2 




Reset numerator & denominator 
accumulators 

NU=DE=0, Reset max/min angles & 
range sum 

Set coordinates of nadir to next 
position 

a = a + (I-l)Aot 
B = + Cl-1)A3 

Y = Yq + (I-I)Ay 

V 

Compute limits of antenna foot- 
print on the ground plane. Reset 
scene line count. 


Set integration limits to cover 
only the antenna foot-print. 
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original page B 

OF POOR QUALITY 



Integration over antenna 
footprint completed 


READ 
UNIT 8 


Read a strip of ground 
scene data 


NR=NR+1 


Increment scene record/ 
line counter 


NR: K1 


Read scene records until 
within antenna footprint 


X=AX(L-1) 


1m=I 1 /I yes 



Y=AY(MM-1) 

+0.5 


Evaluate R,0,<)), 

0..,,,enu,2:de, Rs 







0R/6WV4L 
OF. Poor 


fs 

quality 



Display BT. R, X, Y, XL, XU, YL, YU 
and Satellite Characteristics; 
i»e>, $, 0, B, @pQ> H*o» H, II 




APPENDIX D 
SOFTWARE LISTINGS 


ooo on non or> o rtoonnorsonoo 00000 


ORIGINAL PAGE IS 


STEP 


1 PROGRAM <« 


VEGETATION-CLASS SUB-SCENE TO MODEL CLASS COUNT SUB-SCENE 


INTEGER ASQR 

DIMENSION ITBUFF( 912» 6> , ASQR <6, 6) , CCLAS<6/ 4) , IW0RD<4) 
DIMENSION MSLINE<60e, 2) 

UNIT ASSIGNMENTS FOR THE TI-980: 

B4 == CRT 

B5 = CRT OR CDR (RECORD COUNT INPUT) 

B6 = LP (MSG FILE) 

B7 = INPUT (VEGETATION MAP) TAPE* <1824 BYTE RECORDS 
B8 = OUTPUT (MODEL SCENE) TAPE, 1216 BYTE RECORDS 
B9 = LP OR CRT DIAGNOSTICS PRINT OUT 


RESET CYCLE COUNTER AND END-OF-FILE FLAG 
ICYC=0 
IEND=0 
ICNT=0 
IE0F=0 

OPEN TAPE UNITS FOR INPUT AND OUTPUT 
CALL SVC( 183, 7, 1824, ITBUFF) 

CALL SVC( 184, 7, 1216, MSLINE) 

CALL FOR INPUT DATA 
WRITE (4, 3) 

3 FORMAT (IH , 'INPUT RCD COUNT AND DUMP FLAG < 21 5 ) ' , /, 1 X ) 
READ IN THE NUMBER OF TAPE RECORDS EXPECTED AND DUMP FLAG 

READ(5, 10) IREC, ID 

10 F0RMAT(2I5) 

WRITE(6, 11 ) IREC, ID 

11 FORMATdH , 'EXPECT', 16, ' RODS', 5X, 'DUMP FLAG', 16) 
IRM5=IREC-50 

5 CONTINUE 
ITRCD=0 

IFdEND. LT. 0) GO TO 16 

RESET END OF BUFFER LINES BEFORE READING TAPE 

DO 8 M=l, 6 

DO 7 L=800, 912 

7 ITBUFF(L> M)=32639 

8 CONTINUE 

GET A TAPE RECORD 6-TUPLE, EACH <=912 BYTE PAIRS 
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6 DO 15 N-1.6 

IF( ICNT. EQ, IREC) GO TO 150 
C 

C INPUT A TAPE RECORD FROM UNIT NDR 7 
CALL SVC(183» 1» 1824* ITBUFF< 1, N) ) 

IFdCNT. GT. IRM5) CALL TCK(IEOF) 

IFUEOF. NE. 0) 00 TO 150 
C 

ICNT*ICNT+1 

ITRCD=ITRCD+1 

15 CONTINUE 

16 CONTINUE 
C 

C BUILD 2 LINES OF MODEL SCENE DATA 
DO 125 K=l,304 
I = l-f (K-l)*3 
C 

C BUILD 6-BY-6 WORD-BOX FROM 6-BY-6 VEGETATION BYTE-BOX 
DO 30 U=l,6 
DO 29 L-1, 3 
LK~L~1 
LLsLKX-2-t-l 

TEMP=ITBUFF( I-+LK; J) 

C RIGHT JUSTIFY L-H BYTE 
ASQR ( LL, J ) =TEMP/256 
C FIND VALUE OF R-H BYTE 

ASQR ( LL-H * J ) =TEMP - ( ASQR ( LL , J > »256 ) 

29 CONTINUE 

30 CONTINUE 
C 

C 

C RESET THE COUNTERS FOR CLASSES IN EACH MODEL SCENE BOX 
DO 35 J=l* 4 
DO 34 IK=1, 6 
CCLAS( IK, J)=0 

34 CONTINUE 

35 CONTINUE 
C 

C COUNT CLASSES IN 6-BY-6 BYTE-BOX FOR 4-BY-4 WORD-BOX 
DO 44 JI = 1, 6 
DO 43 KI-1, 6 
IK=Kl/4+l-KJI/4)*2 
C 

C VALUE OF JI INDEXES LINE NBR 
C VALUE OF IK " COLM NOR 

C COUNT NUMBER OF VALUES IN EACH OF 6 CLASSES FOR EACH 
C OF 4 BOXES IN ASQR (INDEXED BY IK) 

CALL COUNT< JI, KI, IK, ASQR, CCLAS) 

43 CONTINUE 

44 CONTINUE 
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BUILD 2 MODEL SCENE WORDS FOR EACH OF TWO LINES 

DO 75 IK*1*4 
IW0RD< IK)»0 

BUILD MODEL SCENE DATA SUB-WORD INDEXED BY IK 
DO 70 KK=1, 5 
KP=KK-1 

IWORD(IK) = lWORD(IK)+IFIX(<CCLAS(KK; IK)/9. )#7. +0. 5)*IFIX<8. *#KP > 
70 CONTINUE 
75 CONTINUE 

PUT 4 SUBWORDS INTO TWO MODEL SCENE LINES 


II = 1 + <K”U*2 
DO 80 LL=1> 2 

MSLINEdI, LL) = IW0RD(2*LL-1 ) 

MSLINEUI-t>l, LL) = IW0RD(2*LL> 

80 CONTINUE 

DIAGNOSTICS *#«■ 

IF(ID. EQ. 0) GO TO 125 

IF< (K. GT. 265), AND. (K. LE. 270). AND. (ICYC. GE. 5). AND. (ICYC. LE. 10) ) 

♦ GO TO 140 

IF< (K. LE. 265). OR, (K. GT. 270). OR. (ICYC. LE. 252). OR. (ICYC. GT. 267) ) 

♦ GO TO 125 

140 CONTINUE 
Ip3=I+2 
IIP1 = I I-fl 

WRITE(‘?* 142) ( (ITBUFF(L. N), L*I, IPO), N=l> 6), ASQR, CCLAS, 

♦ IWORD, ( (MSL1NE(MM, NN). MM=II, I IPl ) , NN=1 , 2) , K, ICYC 
142 FORMATdH , dRCDr /,6(1X,3I7, /)» IX, 'BYTE-BOX', /,6(1X»6I7, /), IX, /, 

$1X, 'CLASS CNT', /, 4(1X, 617, /), IX, /, 

SIX, 'WORD BITS', /, IX, 418, /, 

SIX, 'TAPE RCD WORDS', /, 2( 1 X, 217, / ) , 1 X, /, 

SIX, 'K=', 17, 5X, 'CYCLE*', 17, /, IX) 

*** END OF DIAGNOSTICS 
125 CONTINUE 

ICYC*ICYC+1 

C COMPLETED 2 LINES OF MODEL SCENE 608 WORDS/LINE 
C WRITE 2 LINES TO TAPE UNIT NBR 2 

CALL SVC (184, 3, 1216, MSLINEd, 1) > 

CALL SVC (184, 3, 1216, MSLINEd, 2) ) 

WRITE(4, 141 )ICNT 

141 FORMATdH , 5X, 'RCD CNT', 16) 

IFdEND. LT, 0) GO TO 151 

GO TO 5 

150 IFdTRCD. NE. 0) GO TO 300 
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C COMPLETED TAPE FILE AND MODEL S'JB-SCENE 
151 IRCD*XCYC#6 
ILINESsICYC#2 
ENDFILE 8 

WRITE<6,200) IRCD. ILINES. ICYC, ITRCD, IREC, ICNT 
200 FORMATdOX, . .END OF INPUT TAPE: O 16, ' RECORDS', 

«SX, 16, ' MODEL SCENE LINES', 5X, 16, ' LINE-PAIR CYCLES',/, 
*5X, 16, ' TAPE RECORD RESIDUE',/, 

$5X, 16, ' RECORDS EXPECTED, RECORD C0UNT=',I6) 

STOP 

300 CONTINUE 

FINISH OFF MODEL-SCENE WITH BACKGROUND DATA, >7F7F 
NXL=ITRCD-d 
DO 310 LL«NXL, 6 
DO 309 LC=1, 912 
ITBUFF(LC, LL) =32639 

309 CONTINUE 

310 CONTINUE 
IEND=-1 
GO TO 5 
END 


n > > > > > 


ORIGINAL PAGE « 
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SUBROUTINE TCK(IEOF) 



IDK=0 




DATA 

>C385 

i TURN OFF MPB 

200 

LDM 

?<10 

;OET STATUS WORD 


DATA 

>C38D 



RMO 

M» A 

;MOVE IT TO ACCUM 


STA 

S<23i I 

; STORE AT ISTAT 


CPL 

za 1 

i >8000, COMPLETE? 


SEQ 




BRU 

S<500 

; NOT COMPLETE, CHECK FOR EOF 


BRU 

Zt777 

i OP COMPLETE 

500 

SNZ 

A 

iSKIP IF NOT ZERO STATUS 


BRU 

?<200 

; WAIT FOR CHANGE 


CRA 

8 

iPUT EOF 7-BIT INTO LOC 15 


SEU 

A 

;SKIP IF NOT EOF 


BRU 

S<666 

iEOF FOUND 


CRA 

2 

iNOT EOF, CK PARITY 


SEV 

A 

i SKIP IF NO •• ERROR 


BRU 

S<999 

i SIGNAL ERROR 


BRU 

Z>.777 

; RETURN 

979 

WRITE <4. 

1000) 


1000 

FORMAT (5X* 'TAPE ERROR') 


600 

CONTINUE 




DATA 

>C386 

; TURN ON MPB 


WRITE<6» 

900) I ST AT 


900 

FORMAT (5X 

, 'STATUS', 16) 



BRU 

S<777 

i RETURN 

666 

IE0F=1 



667 

CONTINUE 



DATA 

>C386 

; TURN ON MPB 

800 

RETURN 



777 

IE0F=0 




GO TO 667 


10 

DATA 

:>009C 


11 

DATA 

>8000 


23 

DATA 

ISTAT 



END 
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SUBROUTINE COUNT( J, K, I» ASQ, CCL> 

INTEGER ASQ 

DIMENSION ASQ< 6, 6) ,CCL<6>4)^ 

C *** *** *** *** *** *■«•* * 

c 

C INSERT COUNTING RULES FOR VEGETATION SCENE HERE 
C 

C KERRVILLE KERRVILLE 

C 

IF( <ASQ<K, J). LE. 0). OR. <ASQ<J^, J). GE. 23) ) GO TO 100 
GO TO (1, 2, 3, 4. 5, 6. 7, B. 9, lOi 11. 12, 13. 14. 15, 16, 17. 18, 19. 20. 21. 22), 
«ASQ<K. J) 

100 RETURN 
C 

C WATER — -CLASS 1 

C 

20 CCL< 1. I )=CCL< 1, I )+l 
GO TO 100 

C 

C URBAN S< BARE SOIL -CLASS 2. 3 

C 

5 CONTINUE 

21 CONTINUE 

22 CCL(2. I )=CCL(2. I ) + l 
GO TO 100 

C 

C MIXED SOIL AND VEGETATION CLASS 4 

C 

17 CONTINUE 

18 CONTINUE 

19 CCL<4, I >=CCL(4. I )+l 
GO TO 100 

C 

C FULLY VEGETATED (NON-FOREST) CLASS 5 

C 

1 CONTINUE 

2 CONTINUE 

3 CONTINUE 

4 CONTINUE 

6 CONTINUE 

7 CCL(5. I )=CCL(5, I )+l 
GO TO 100 

C 

C FOREST LAND CLASS 6 

C 

8 CONTINUE 

9 CONTINUE 

10 CONTINUE 

11 CONTINUE 

12 CONTINUE 
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C *** AUSTIN AUSTIN *#* 

C 

IF( <ASQ(K» J). LE, 0), OR. (ASQ(K, J). GE. 34) ) GO TO 100 
GO TO < 1» 2, 3. 4, 5, 6, 7i 8i 9, 10, 1 li 12, 13, 14/ 15, 16» 17, 18, 19i 20» 
*21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33), ASO(K, J) 

100 RETURN 
C 

C WATER CLASS 1 

C 

25 CONTINUE 

26 CONTINUE 

27 CCL( 1, I )=CC.La, I ) + l 
GO TO 100 

C 

C URBAN ?< BARE SOIL CLASS 2, 3 

C 

19 CONTINUE 

28 CONTINUE 

29 CONTINUE 

30 CONTINUE 

31 CONTINUE 

32 CCL(2, I )=CCL(2, I )+l 
GO TO 100 

C 

C MIXED SOIL VEGETATION CLASS 4 

C 

10 CONTINUE 

11 CONTINUE 

15 CONTINUE 

16 CONTINUE 

17 CONTINUE 

18 CONTINUE 

20 CONTINUE 

21 CCL(4, I )=CCL(4, I ) + l 
GO TO 100 

C 

C FULLY VEGETATED (NON-FORESTED) CLASS 5 

C 

1 CONTINUE 

2 CONTINUE 

3 CONTINUE 

5 CONTINUE 

6 CONTINUE 

7 CONTINUE 

8 CCL<5, I )=CCL(5, I)+l 
GO TO 100 

C 

C FOREST LAND CLASS 6 

C 

4 CONTINUE 
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9 CONTINUE 

12 CONTINUE 

13 CONTINUE 

14 CONTINUE 
22 CONTINUE 

33 CONTINUE 

34 CCL<6, I )=CCL(6, I )+l 

3 CONTINUE 

4 CONTINUE 
GO TO 100 

C 
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*** BROWNWODD DRQWNWOOD 

IF( (ASQ<K* J). LE. 0). OR. <ASQ<K» J). GE. 8) ) GO TO 100 
GOTO < 1, 2, 3, 4, 5, 6, 7), ASQ<Ki J> 


100 RETURN 

WATER CLASS 1 

6 CCL( 1, I )=CCL( 1/ 1)4-1 
GO TO 100 

URBAN ?< BARE SOIL CLASS 2, 3 

7 CCL(2, I )=CCL<2. I )+l 
GO TO 100 

MIXED SOIL !< VEGETATION CLASS 4 

5 CCL(4, I)-CCL(4, I)+l 
GO TO 100 

FULLY VEGETATED (NON-FOREST) CLASS 5 


1 CONTINUE 

2 CCL(5, I )-CCL(5, I )+l 
GO TO 100 

FOREST-LAND CLASS 6 

3 CONTINUE 

4 CCL(6, I )=CCL(6, I ) + l 
GO TO 100 

END 
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«#«• WACO WACO *** 

IF< (ASQ<K, J). LE, 0). OR. (ASQ<K, J). GE, 9) ) 00 TO 100 


GO TO ( I, 2, 3» 4» 5» 6. 7, 8>, ASQ(K, J) 

100 RETURN 

WATER -CLASS 1 

7 CCL< 1, I )=CCL< 1, I )+l 
GO TO 100 

URBAN !< BARE SOIL CLASS 2i 3 

8 CCL(2, I )=CCL<2, I)+l 
GO TO 100 

MIXED SOIL ?< VEGETATION -CLASS 4 

6 CCL<4, I )=CCL(4. D-^-l 
GO TO 100 

FULLY VEGETATED (NON-FOREST) CLASS 5 


1 CONTINUE 

2 CCL(5. I)=CCL(5, I )+l 
GO TO 100 

FOREST-LAND CLASS 6 

3 CONTINUE 

4 CONTINUE 

5 CCL(6i I )=CCL(6, I)+l 
GO TO 100 


END 
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*** BRYAN BRYAN *** 


IF( <ASQ(K. J). LE. 0). OR. (ASQ(K, J). GT, 39) )G0 TO 100 

GO TO <1.2. 3. 4* 6, 7. Q, 9, 10» 11« 12^ 13. 14. 15. 16. 17. 18. 19. 20. 

#21. 22. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. 33. 34. 35. 36. 37. 38. 39). ASQ<K. J) 
100 RETURN 

OPEN WATER CLASS 1 

24 CONTINUE 

26 CCL( 1. I >=CCL< 1. I >+l 
GO TO 100 

URBAN AND DARE SOIL CLASS 2. 3 

8 CONTINUE 

21 CONTINUE 

22 CONTINUE 

23 CONTINUE - 

33 CONTINUE 

34 CONTINUE 

38 CCL<2. I >=CCL<2, I)+l 
GO TO 100 

MIXED SOIL AND VEGETATION 

MIXED SOIL AND VEGETATION ‘ - - -CLASS 4 

9 CONTINUE 
13 CONTINUE 
15 CONTINUE 

18 CONTINUE 

19 CONTINUE 

27 CONTINUE 

29 CONTINUE 

30 CONTINUE 

31 CONTINUE 

35 CONTINUE 

36 CONTINUE 

37 CONTINUE 

39 CCL(4. I)=CCL(4. I)+l 
GO TO 100 

FULLY VEGETATED < NON-FOREST) CLASS 5 

1 CONTINUE 

3 CONTINUE 

4 CONTINUE 

5 CONTINUE 
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6 CONTINUE 
14 CONTINUE 

32 CCL<5, I )=CCL<5, I )+l 
GO TO 100 

FOREST-LAND CLASS 6 

2 CONTINUE 

7 CONTINUE 

10 CONTINUE 

11 CONTINUE 

12 CONTINUE 

16 CONTINUE 

17 CONTINUE 
20 CONTINUE 
25 CONTINUE 

28 CCL<6, I >=CCL<6. I)+l 
GO TO 100 

END OF BRYAN RULES 
END 
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C CORSICANA CORSICANA *** 

C 

IF( (ASQ(K, J). LE. 0). OR. <ASQ(K» J>. OT. 36) )GO TO 100 
GO T0( 1* 2. 3, 4. 5. 6 , 7 , 8, 9, 10, 11, 12, 13. 14, 15, 16, 17, 18, 19, 20, 
$21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36), ASQ(K, J) 
100 RETURN 
C 

C WATER CLASS 1 

C 

17 CONTINUE 

18 CONTINUE 

19 CONTINUE 

20 CCL( 1, I )=CCL< 1, I ) + l 
GO TO 100 

C 

C URBAN BARE SOIL CLASS 2, 3 

C 

29 CONTINUE 

30 CONTINUE 

31 CONTINUE 

32 CONTINUE 

34 CONTINUE 

33 CONTINUE 

35 CONTINUE 

36 CCL(2, I )=CCL(2, I ) + l 
GO TO 100 

C 

C MIXED SOIL >!•< VEGETATION CLASS 4 

C 

21 CONTINUE 

22 CONTINUE 

23 CONTINUE 

24 CONTINUE 

25 CONTINUE 

26 CONTINUE 

27 CONTINUE 

28 CCL(4, I )=CCL(4, I)+l 
GO TO 100 

C 

C FULLY VEGETATED (NON-FOREST) — —CLASS 5 

C 

1 CONTINUE 

2 CONTINUE 

3 CCL<5, I)=CCL<5, 1) + 1 
GO TO 100 

C 

C FOREST LAND CLASS 6 

C 

4 CONTINUE 

5 CONTINUE 
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6 CONTINUE 

7 CONTINUE 

8 CONTINUE 

9 CONTINUE 

10 CONTINUE 

11 CONTINUE 
IS CONTINUE 

13 CONTINUE 

14 CONTINUE 

15 CONTINUE 

16 CCL(6, I )*CCL<6, I )4.1 
GO TO iOO 

END OF RULES FOR CORSICANA 
END 
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C *** HOUSTON HOUSTON 
C 

IF( <ASQ(K, U). LE. 0). OR. (ASQ<K, J). GT. 39) )GO TO 100 
GO TO (1, 2» 3, 4, 5, 6, 7, 8. 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 
i23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39), ASQ<H, 
100 RETURN 
C 

C WATER CLASS 1 

C 


34 CCLd, I )=CCL< 1, I) + l 
GO TO 100 

35 GO TO 34 
37 GO TO 34 

C 

C *** URBAN BARE SOIL 
C 

25 CCL(2, I >=CCL(2, I)+l 
GO TO 100 

26 CONTINUE 

27 CONTINUE 

28 CONTINUE 

29 CONTINUE 

30 CONTINUE 

31 CONTINUE 

32 CONTINUE 

33 CONTINUE 
39 GO TO 25 

C 

G MIXED SOIL VEGETATION 

15 CONTINUE 

16 CONTINUE 

17 CONTINUE 

18 CONTINUE 

19 CONTINUE 

20 CONTINUE 

21 CCL<4, I >=CCL(4, I )'<-l 

22 GO TO 100 
C 

C FULLY VEGETATED (NON/FOREST) 
C 

1 CONTINUE 

2 CONTINUE 

3 CONTINUE 

23 CONTINUE 

24 CCL(5, I )=CGL<5, I)+l 
GO TO 100 

C 

C FOREST-LAND 

4 CONTINUE 

5 CONTINUE 


CLASS 2, 3 


-CLASS 4 


CLASS 5 


CLASS6 
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6 CONTINUE 

7 CONTINUE 

8 CONTINUE 

9 CONTINUE 
ID CONTINUE 

11 CONTINUE 

12 CONTINUE 

13 CONTINUE 

14 CONTINUE 
36 CONTINUE 

38 CCL<6, I )=CCL<6i I )+l 
GO TO 100 

END 
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! 


c 

c 

c 



c 

c 

c 


c 

c 

c 


c 

c 


c 

c 

c 


c 

c 

c 


c 


♦#* LUFKIN LUFKIN ### 


IF( <ASQ(Ki U>. LE. 0). OR. <ASQ(K> J). 5T. 14) ) GO TO 100 
GO TO < 1* 2, 3, 4. 5. 6, 8< 9, 10/ 111 12, 13, 14), ASQ(Ki J) 

100 RETURN 

WATER - -CLASS 1 

10 CCL(1, I )=CCL(1, I)+l 
GO TO 100 


URBAN t< BARE SOIL -CLASS 2» 3 

11 CCL<2i I )=CCL(2/ I )‘i-l 

C CLASS 3 = SAME AS CLASS 2 
GO TO 100 

MIXED SOIL VEGETATION ---------- -CLASS 4 

12 CCL(4, I )=CCL(4, I )+l 
GO TO 100 

13 GO TO 12 

14 GO TO 12 

VEGETATED < NON-FOREST) CLASS 5 

1 CCL<5, I )=CCL<5, I )+l 
GO TO 100 

2 GO TO 1 

FOREST- CLASS 6 

3 CCL<6/ I )-CCL<6, I )+l 
GO TO 100 

4 GO TO 3 

5 GO TO 3 

6 GO TO 3 

7 GO TO 3 

8 GO TO 3 

9 GO TO 3 
END 


r 
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FILL TAPE DATA SCENE WITH BACKGROUND 


ASSIGN B7 TO MT 
ASSIGN B4 CRT 

INTEGER I SCENE (2496) 

C FILL SUPER-SCENE LINE BUFFER WITH BACKGROUND 
DO 12 I C= 1,2496 
ISCENE(IC) =0 
12 CONTINUE 

C OPEN TAPE UNIT FOR 2496 WORDS PER RECORD 
CALL SVC (183, 7, 4992, ISCENE( 1 ) ) 

ICNT=0 

WRITE(4, 10) ICNT 

10 FPRMATUH , 'TAPE UNIT OPEN, RCD CNT=',I6) 

DO 24 KR=1, 1650 

C OUTPUT A RECORD FROM THE BUFFER FILLED WITH ZERO 
CALL SVC( 183, 3, 4992, ISCENE( 1 ) ) 

ICNT=ICNT+1 
24 CONTINUE 

C OUTPUT TWO END OF FILE MARKS ON THE TAPE 
CALL SVC( 183, 10, 4992, ISCENE( 1 ) ) 

WRITE(4, 30) ICNT 

30 FORMATdH , 'TAPE UNIT CLOSED, RCD CNT=',I6) 
STOP 
END 
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«#«##«■*«-# STEP 2 PROGRAM ##*###### 

MODEL SUPERSCENE FROM MODEL SUBSCENES 

UNIT ASSIGNMENTS: 

8 = INPUT MODEL SCENE 

9 = INPUT SUPERSCENE 
10 = OUTPUT SUPERSCENE 

ISCN = (INPUT SCENE SEQ NBR: 1 THROUGH B> 

1 = LUFKIN 

2 = HOUSTON 

3 = CORSICANA 

4 = BRYAN 

5 = WACO 

6 = AUSTIN 

7 - BROWNWOOD 

8 = KERRVILLE 

IMSL = NBR OF LINES IN THE SCENE BEING PROCESSED 
SYSTEM OPERATION: 

1. READS UNIT 5 FOR OPERATION PARAMETERS 

2. OVERLAYS SUPER“SCENE WITH MODEL SUB -SCENE SPECIFIED BY ISCN 
THREE TAPES SHOULD BE MOUNTED, INPUT = MODEL SUBSCENE 

INPUT = OLD SUPERSCENE 
OUTPUT = NEW SUPERSCENE 

INTEGER#2 ISCENE, MSCENE 
DIMENSION ISCENE(2496, 2 > , MSCENE ( 608, 2) 

COMMON/FILLER/ ISCENE 

GET RUN PARAMETERS FROM CARD READER 

READ (5, 10) ISCN, IMSL 

10 F0RMAT(3I5) 

BACKGROUND FILLED WITH >0000 
WRITE<6, 12) ISCN, IMSL 

12 FORMAT( IHl. lOX, 'STEP 2 PROGRAM RUN PARAMETERS: 

#20X, 'INPUT SCENE=', 16, lOX, 'INPUT SCENE LINE COUNT= ' , 16, / ) 

CHECK MODEL SUB-SCENE SEQ NBR 

15 CONTINUE 

GOTO (100,200,300,400,500,600,700,800) , ISCN 

THIS SECTION FOR LOADING MODEL SUBSCENE 
100 CONTINUE 
C 

LC0L=188B 
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LNBR«0 

LN=0 

WRITE<6, 110) ISCN» LNBR» LCOL 

110 FQRMAKIH i 10X» 'BEGINNING SCENE', 16, LUFKIN AT LINE', 16, 
COLUMN', 16, /) 

C 

225 CONTINUE 

C GET LINE FROM MDDEL-SCENE INTO BUFFER PAIR NBR 1 
READ (8, 226, E<ND=1777) <MSCENE(M, 1 ), M=l, 608) 

226 FORMATS 19 (32A2) ) 

GET LINE FROM SUPER-SCENE INTO NBR 1 BUFFER PAIR 
READ<9, 227M ISCENEdl, 1), 11 = 1, 2496) 

227 FORMAT < 39 < 64A2 ) ) 

OVER-LAY 1ST SUB-SCENE LINE ONTO 1ST SUPER-SCENE LINE 

DO 230 MC= 1,608 
IC=LC0L+MC-1 

IF(1SCENE< IC, 1). EQ. 0)ISCENE(IC, 1 )=MSCENE(MC, 1) 

230 CONTINUE 

OUTPUT THE REFORMATTED SUPER-SCENE LINE FROM BUFFERS 1 

WRITE( 10, 227) (ISCENEdl, 1), 11 = 1,2496) 

LNBR=LNBR+1 
LN=LN+1 

IF(LN .GE. IMSL) GO TO 1777 
GET ANOTHER PAIR OF LINES INTO NBR 2 BUFFERS 

READ (8, 226, END=1777 ) ( MSCENE ( M, 2), M=l, 608) 

READ(9, 227) < ISCENEd I, 2), 11 = 1, 2496) 

OVER-LAY 2ND SUB-SCENE LINE ONTO 2ND SUPER-SCENE LINE 
DO 240 MC=1,608 
IC=LC0L+MC-1 

IFdSCENEdC, 2). EQ. 0) ISCENE( IC, 2) -MSCENE <MC, 2) 

240 CONTINUE 

C OUTPUT THE REFORMATTED SUPER-SCENE LINE FROM BUFFERS 2 
WRITEdO, 227) (ISCENEdl, 2), 11 = 1, 2496) 

LNBR=LNBR+1 

LN=LN+1 

C COMPLETED TWO LINES INTO THE NEW SUPER-SCENE 
IF(LN .GE. IMSL) GO TO 1777 
GO TO 225 
C 

C FINISH WRITING THE SUPER-SCENE ONTO NEW TAPE 
C POINT TO THE NEXT LINE TO PROCESS 
1777 LAST = LNBR . 

1800 IF(LNBR. GE. 1650) GO TO 1801 

READ (9, 227) (ISCENEdl, 1), 11 = 1,2496) 
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WRITEnO,227) (ISCENE<II, 1), IIal>2496> 

LNBR«LNBR-fl 
GO TO 1800 
1801 CariTINUE 

SUPER-SCENE NOW ON (NEW) TAPE 
ENDFILE 10 

WRITE(6, 1900) LAST, LNBR,LN 

1900 FORMATdOX, '—LAST MODEL SCENE LINE=',I6, ' ON SUPERSCENE', 

$/, MX, ' LAST LINE NBR-', 16, ' LAST MODEL SCENE LINE=', 16, /) 

STOP 

FOR LOADING HOUSTON SUBSCENE ONTO SUPER SCENE 
200 LNBR=680 
LC0L=1848 

WRITE<6, 210) ISCN, LNDR, LCOL 

210 FORMATdH , lOX, 'BEGINNING SCENE', 16, ', HOUSTON AT LINE', 16, 
COLUMN', 16, /) 

250 CALL SKIP(LNBR) 

INITIALIZE MODEL SCENE LINE NUMBER 
LN=0 

GO TO 225 

FOR LOADING SUB-SCENE FROM CORSICANA 
300 LNBR=99 
LC0L=1330 

WRITE(6,310) ISCN, LNBR, LCOL 

310 FORMATdH , lOX, 'BEGINNING SCENE', 16, ', CORSIC. AT LINE', 16, 
COLUMN', 16, /) 

GO TO 250 

FOR LOADING SUB-SCENE FROM BRYAN 
400 LNBR=780 
LC0L=1292 

WRITE<6, 410) ISCN, LNBR, LCOL 

410 FORMATdH , lOX, 'BEGINNING SCENE', 16, ', BRYAN AT LINE', 16, 
COLUMN', 16, /) 

GO TO 250 

FOR LOADING SUB-SCENE FROM WACO 
500 LNBR=212 
LC0L=772 

WRITE<6, 510) ISCN, LNDR, LCOL 

510 FORMATdH , lOX, 'BEGINNING SCENE', 16, ', WACO AT LINE', 16. 

*', COLUMN', 16, /) 

GO TO 250 

FOR LOADING SUB-SCENE FROM AUSTIN 
600 LNDR=892 
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LC0L=734 

WR1TE(6>610) 1SCN» LNDR. LCOL 

610 FORMATdH , 10 %, 'BEGINNING SCENES 16. ', AUSTIN AT LINES 16. 
#S COLUMNS 16. /> 

GO TO 250 

FOR LOADING SUB-SCENE FROM BROWNWOOD 
700 LNBR=281 
LC0L=216 

WRITE<6. 710) ISCN. LNBR, LCOL 

710 FORMAT <1H . lOX. 'BEGINNING SCENE'. 16. ', BROWNW. AT LINE'. 16. 
COLUMN'. 16. /) 

GO TO 250 

FOR LOADING SUB-SCENE FROM KERRUILLE 
800 LNDR=963 
LC0L=170 

WRITE(6.810) ISCN. LNBR. LCOL 

010 FORMATdH . lOX. 'BEGINNING SCENE'. 16. '. ;^ERRUI. AT LINE'. 16. 
COLUMN', 16, /) 

GO TO 250 
C 

END 

SUBROUTINE SKi'P(LCNT) 

INTEGER*2 ISCENE<2496. 2) 

COMMON/FI LLER/ISCENE 
C 

DO 310 K--=1,LCNT 

READ (9, 10) ( ISCENE( J, 1 ), J=l, 2496) 

WRITEC 10. lOX ISCENECL. 1 ). L=--l. 2496) 

10 FORMAT ( 32 ( 78A2 ) ) 

310 CONTINUE 
RETURN 
END 
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KJDEL OF ORBITING MICROWAVE RADIOMETER 

I'^PUT 3 CARDS FOR EACH PASS; 

CARD #1: 

SF5. 0 — AZIM, BMWID, INC ID* ALT. TEMP 
DEG DEG DEG KM DEG C 
215 — DOWN ?< CROSS RANGE SCENE SIZE 
CARD#2; 

4F5. 0— SOIL MOIST X Y DELTA 
ROUGHNESS X ?< Y DELTA 
215 — NADIR MOTION STEP SIZE. 

BAND IDENTIFICATION # 

2F6. 1— SOIL MOIST ROUGHNESS 

CARD #3; 

2F5. 0 — X !< Y STARTING COORD. O'" ANTENNA FOOT PRINT 


INPUT UNITS: 

5=CDR 

7=MT1 (MOIST, ROUGH) 
8=MT2(VEG CLASS) 

OUTPUT UNITS: 

fe=LP(BT‘ HISTORY) 

9sMT3(BT HISTORY FARMS) 


C 


SYSTEM BAND: 1=L. 2= 

AZIM-AZIMUTH(=>0, =< 
BWID-BEAMWIDTH 
ANGI-INCIDENT ANGLE 
(=>0, =<45 ) 

TARGET COORDINATES: 
TOP-X COORDINATE 
CENT-Y COORDINATE 
■SFC-SURFACE 


C. 3=X 

90) ALT-SATELLITE ALTITUDE 

TEMP-BACKGROUND TEMP. (DEG C) 
SCENE BOUNDARY LIMITS: 
LMDR-DOWNRANGE 
LMCR-CROSSRANGE 
SCENE INTEGRATION STEPS: 
DOWN-DOWNRANGE STEP 
CROS-CROSSRANGE STEP 
ELEV-ELEVATION STEP 


NADIR DISPLACEMENT 
DNAD-DOWNRANGE 
CNAD-CROSSRANGE 
ENAD-ELEVATION 


DOU^'LE PRECISION RSUM, RCNT 
INTEl*^R*2 ILINE, IMOIST 

DIMENS. HN ILINE(2496), IMOIST ( 2496 ) . B Y < 325 ) . BR ITV<325 ) . BR ITH( 325 ) 
EQUIVALEn"E(X, TOP). (Y, CENT). (ANGI. ROLL) 

COMMON/COORu.'V, XL. XU. Y. YL, YU, XI. X2. Yl, Y2 
COMMON/DAT/ILINt^. ''HOIST 
COMMON/ROUGH/RH. RV 
COMMON/GAN/BWID, XPI. IF 
COMMON/OR lEN/ROLL, AZIM. ALT, A. B 
COMMON/TRIG/COSR. SINR. COSA, SINA 
COMMON/CLAS/CLASUM ( 6 ) 

CONRAD ( ANGLE ) =ANGLE/57. 2957 
IF=3 

XPI=3. 14159 
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C INITIALIZE RECORD NBR 

C SCAN PARAMETERS READ IN» VARIABLES INITIALIZED 
IRCD*0 

200 CALL INPUT<AZIM.BUIID,ANGl»LMDR,LMCR, ALTO, TEMP) 

READ MOISTURE ROUGHNESS RATES, INTEGRATION STEP SIZE. BAND 
READ < 5, 8)DSMY, DSMX, RUPY, RUFX, I STEP, NBND, SO ILK, ROUK 

INITIALIZE SATELLITE NADIR COORDINATES 
READ ( 5, 1 3 ) X, Y, RESL, DNAD 
13 FORMAT <3F5. 0, F5. 4 ) 

INITIALIZE SURFACE ELEVATION(MSL) 

SFC0=0. 0 

SET DOWN AND CROSS RANGE INTEGRATION UNIT SIZE, SET ELEV. UNIT STEP 
D0WN«1. 0 
CR0S=1. 0 
ELEV-0. 0 

INITIALIZE ALTITUDE FOR SPECIFIED RESOLUTION 
IF(RESL. NE. 0. )CALL CALT( RESL, ALTO) 

SET NADIR DISPLACEMENT UNIT SIZES 
CNAD-l. 0 
ENAD^O. 0 
DELZ-0. 0 

INITIALIZE LINE COUNTER, 

AND SET DELTA-THETA STEP 
11=0 

INITIALIZING COMPUTATIONAL VARIABLES 
DA=DOWN*CROS 

SET UNIT SCALE FACTOR SCALE ALTITUDE 
S=0, 24 

salt=alto 

ALTO=ALTO/S 

WRITE<6, 2)AZIM, TEMP, BWID, ANGI, LMDR, LMCR. SALT, TOP, DOWN, 

$DNAD, CENT, CROS, CNAD, SFCO, ELEV, ENAD 
WRITE<9, 210)AZIM, TEMP, BWID, ANGI, LMDR, LMCR, SALT, TOP, DOWN, 

♦DNAD, CENT, CROS, CNAD, SFCO, ELEV, ENAD 
ANGSV=ANGI 
C CONVERT ANGLES TO RADIANS 
AZSV=AZIM 
BWSV=BWID 
TPSZ=TEMP 
ANGI =CONRAD( ANGI ) 

AZIM=CONRAD(AZIM) 

BWID«CONRAD<DWID) 

COSR*COS(ROLL) 

S1NR»S1N(R0LL) 

COSA*COS(AZIM) 
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C EVALUATE SURFACE RANGE* AND INITIALIZE NADIR COORDINATES 
SFRG=ALTO*TAN ( ROLL ) 

SRSINA«SFRC*SINA 

AO«TOP-fSRSINA 

SRCOSA«SFRG#COSA 

BO«CENT-SRCOSA 

IF( ISTEP, NE. 0)G0 TO 53 

ISTEP=IFIX(ALT0*TAN<BWID/2. )) 

53 WRITE <6* 3>SFRG. AO* DO* SO ILK* ROUK 
WRITE (6* 9) ISTEP* NBND, RESL 

C BEGIN COMPUTATION #'*##**##*#>##***********#-»*i^-!<^**<»****-»#<HHHH»**<HHHH» 
DO 10 J=l, LMDR. ISTEP 

C RESET NUMERATOR h DENOMINATOR ACCUMULATOR 
D*sO. 0 
XNH-0. 0 
XNV=0. 0 
11 = 11+1 
DO 57 LL=1< 6 
57 CLASUM(LL)=0. 0 

C INCREMENT NADIR COORDINATES TO NEXT POSITION 
L=J-1 

A=AO+L*DNAD 
B=B0+L»CNAD ’ 

C PUT NADIR POINT IN PLOT BUFFER 
BY( 1 1 )=B+SRCOSA 

C STEP ALTITUDE AND SURFACE 
ALT=ALTO+L*ENAD 
SFG=SFCO+L*DELZ 

C INITIALIZE RANGE SUMMATION REGISTER 
RSUM=0. 0 

C FIND UPPER AND LOWER LIMITS OF ANTENNA SPOT SIZE 
X=A-SRSINA 
Y=BY(II > 

CALL XYLM1T(DX. DY) 

IF( J- NE, DGO TO 101 
WRITE(6, 12)XU. YU* XL* YL, DX. DY 

WRITE<9, 211 )SFRGi AO, BO, SOILK, ROUK, ISTEP* NBND* RESL, 
eXU, YU, XL, YL, DX, DY 
WRITE<6, 5) 

101 CONTINUE 

C SET LIMITS ON INTEGRATION IN X-DIRECTION 
K1 = IFIX<XL-. 5) 

IF(K1. LE. 0)K1 = 1 
IF(K1.GT. 1650)Kl = 1650 
K2=IFIX(XU+. 5) 

IF(K2. LT. 1)K2=1 
IF<K2. GT. 1650)K2*1650 

C SET LIMITS ON Y-DIRECTION 
I1*IFIX(YL-. 5) 
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IF( II, LE. 0)11-1 

IF< II. OT. 2476)11*2496 

I2=IFIX(YU+. 5) 

IF(I2. LE. 0)12=1 
IF(I2. OT. LMDR)00 TO 110 
C INITIALIZE RANGE AVERAGING COUNTER 
C SKIP DOWN GRND SCENE TO FOOT-PRINT 
RCNT«0, 0 

IFdRCD.LT. <K1-D) CALL ISKIP ( IRCD, K1 ) 

IFdRCD. GT. (K1'-1))CALL IBACKdRCD. K1 ) 

C SCAN DOWN FDQT-PRINT. LINE-BY-LINE 
DO 30 K=K1,K2 

INPUT SCENE DATA 

READ STRIP OF DATA FROM GROUND SCENE AND MOISTURE OVERLAY 
READ (7/ 300) IMOIST 
READ(8» 300, END=1 10) ILINE 
C INCREMENT OVER ANTENNA-FOOT-PRINT 
IRCD=IRCD-»-l 
T0P=<K-1 )*D0WN+0. 5 
C SCAN ACROSS FOOT-PRINT 
DO 20 1=11, 12 

C SCAN SCENE, ACCUMULATE BRIGHTNESS COMPONENTS 
CENT=d-l )*CROS+Q, 5 
C COMPUTE RANGE TO GROUND CELL UNIT 

RVEC ^SQRT ( ( TOP-A ) #*2+ ( CENT-B ) **2+ ( SFC-ALT ) **2 ) 

C EVALUATE THETA, ZEN, ANT. ANGLE 

T=ARCOS( ( ( TOP-A )*S I NR»SlNA-( CENT-B)* 

«S I NR*COSA-^ < SFC-ALT ) *COSR > /RVEC ) 

C EVALUATE PHE, ANTENNA AZIM. ANGLE 
XDP=(X-A)*COSA-KY-B)*SINA 
YDP= ( Y-B ) *COSA- ( X-A ) *S INA 
P=ATAN<XDP/YDP) 

C COMPUTE ACTUAL INCIDENCE ANGLE TO GND CELL UNIT 
TN=ARCOS(ALT/RVEC > 

C " UNIT AREA, RANGE SUM, RANGE COUNT 
RSUM=RSUM-»-RVEC 
RCNT=RCNT+1. 0 

C EVALUATE ANTENNA TEMPERATURES 

SM=B0ILK*(1, 0+DSMY*Y/2496. +DSMX*X/ 1650. ) 

ROU=ROUK*d. 0*RUFY*Y/2476. +RUFX*X/ 1650. ) 

CALL BR1GHT<SM, ILINE d ) . NBND, TEMP, ROU, BTVSO, BTHSO) 
CALL BCORR < BTVSO, BTHSO, TN, BTV, BTH ) 

TBV=BTV 

TBH=BTH 

BTH=BRTMP ( TBH, TB V, P ) 

BTV=BRTMP (TBV, TBH, P ) 

SAVE=G ( T, P ) *COS < TN ) *DA/RVEC**2 

XNV=XNV+BTV*SAVE 

XNH=XNH+BTH*SAVE 
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D«D+SAVE 
20 CONTINUE 

30 CONTINUE \ 

RAVO»SNGL(RSUM/RCNT)*S j 

DO 36 LL^1< 6 I 

36 CLASUri(LU“-CLASUI1(LL>/RCNT f 

C COMPUTE AND PRINT BRIGHTNESS VALUES 
BRITV( 1 1 >*XNV/D 
BRITH<II>»XNH/D 

WR1TE(6, 4) J, SALT, RAVG, BRITVC 1 1 >, BRITH( 1 1 ) , A, BY< 1 1 ) , CLASUM 
WRITE (9, 400 ) J, SALT. RAVG, BRITV< 1 1 ) , BRITH< 1 1 > , A, BY< 1 1 ) , CLASUM 
C 

10 CONTINUE 
110 WRITE (6, 7) 

ENDFILE9 

GO TO 200 L 

2 FORMATS 'I', ^SATELLITE ATTITUDE «< CHARACTERISTICS'/' AZIMUTH<HEAC 

DING) 3X, F4. 0< 'DEG',15X, 'BACKGROUND TEMP'.FS, 0, 'DEG'/' BEAMWIDTH 
U',10X,F4.0, 'DEG',16X, 'SCENE BOUNDARY LIMITS: '/' INCIDENT ANGLE', I 
H5X, F4. 0, 'DEG', 16X, 'DOWNRANGE', 7X, 14/ ' ',9X, ^ 

@20X, 'CROSSRANGE', 6X, 14/ ' 27X, 

, 'SATELLITE ALTITUDE ',F4. 0, 'KILOMETERS',/, 

IIX,/' TARGET COORDINATES', IIX, 'SCENE INTEGRATION STEPS', lOX, 'NADIR 
R DISPLACEMENT STEPS'/' INITIAL X ',F7. 1, SX, 'DOWNRANGE STEP', 

, 4X, F3. 1, 12X, 'DOWNRANGE STEP ',F5.3/' INITIAL Y' 

,,4X,F7. 1, 8X, 'CROSSRANGE STEP ', F3. 1 , 12X, 'CROSSRANGE STEP ', 

, F3. 1/ ' INITIAL Z ', 4X, 

,F7.1, 8X, 'ELEVATION STEP ', F3. 1 , 12X, 'ELEVATION STEP ',F3, 1) 

3 FORMAT ( 'OSURFACE RANGE ',F6. 1,5X, ' NADIR COORDINATES <',F7. 1, ', ', 
eF7. 1, ') ', 5X, 'SOIL MOIST, CONST. ',F6. 1, 5X, 'ROUGHNESS CONST, ',F6. 2> ] 

4 FORMAT( IH , 15, 7X, F4, 0, 8X, F9, 2, lOX, F8, 2, 4X, F8. 2, 7X,2F8. 2, IX. 

e6(F6. 2, IX) ) I 

5 FORMAT('ORNG STEP 2X, 'ALTITUDE(KM) ', 4X. 'RANGE(KM) ', 6X, 1 

e 'BRIGHTNESS VER/HOR (DECK)', 7X, 'BEAM COORD( X, Y> ', lOX, S 

e'AVG CLASS PERCENTS') 1 

7 FORMAT (IH , 'END OF SCENE',//) | 

8 FORMAT (4F5, 0, 215, 2F6. 1 ) 

9 FORMAT! 'OSTEP SI ZE- ' , 12, 5X, 'BAND IDENT. ', 13, 5X, 'RESOLUTION= ' , F5. 1 )l 

12 FORMAT! 'O', 'XU=', F9. 2, 2X, 'YU=', F9. 2, 2X, 'XL='. F9. 2, 2X, 'YL*', F9. 2, | 

,2X, 'ANTENNA SPOT 5TZE',F7. 2, ' BY',F7. 2) 

210 FQRMAT!4F4. 0, 214, F4. 0, 3!F7. 1, 2F3. 1 ). 65X) 

211 F0RMAT!4F6. 1, F6. 2, 213, F5. 1, 6F7. 2, 49X) ‘ 

300 FORMAT ! 64 ! 39A2 ) ) ! 

400 F0RMAT!I5, 6!FB. 2, 2X), 6F6. 2, 3iX) 1 

END I 

SUBROUTINE INPUT ( A, B, I , D, C, H, T ) I 

REAL I 
INTEGER D, C 

READ(5, 10. ENDs20)A, B, I,H, T, D. C 
10 FORMAT ( 5F5. 0, 215) 
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RETURN 
20 ENDFILE9 
STOP 
END 

FUNCTION TAN(X> 

TAN«SIN<X)/COS(X) 

RETURN 

END 

FUNCTION G(T<P) 

COMMON/GAN/BWID, XPI» IF 
X=( (XPI-T>*XPI >/(BWID/2. ) 

IF< X,EQ. 0. ) GO TO 7 
G*ABS( (SIN(X)/X)**IF) 

5 RETURN 
7 G*l. 0 
GO TO 5 
END 

SUBROUTINE rSKIPUR/K,' 

INTEGER*2 ILINE(a-:^96) . IMaiST<2496) 

COMMON/DAT/ ILINEi I MO I ST 

KM“K-i 

5 IF( IR. EQ. KM) GO TO 7 
C READ(7, 20)IM0IST 

READ (8, 20.END=30)ILINE 
20 FORMAT ( 64 < 39A2 ) ) 

IR=IR+1 
GO TO 5 
7 RETURN 

30 WRITE<6, 35) IR, K 
35 FORMAT ( ' IR= S 16, lOX, 'K= ^ 16) 
ENDFILE9 
STOP 
END 

SUBROUTINE IBACK(IR,K) 

IN7EGER*2 ILINE(2496), IM0IST(2496) 

COMMON/DAT/ ILINE, IMOIST 

KM=K-1 

5 IF( IR. EQ. KM)GO TO 7 
C BACKSPACE 7 

BACKSPACE 8 
IR=IR-1 
GO TO 5 
7 RETURN 
END 
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SUBROUTINE BRIGHKSM# ICW, ND. TP» ROUi BTVSOi DTH50) 

INTEGER#2 I CM 
COMMON/ROUOH/RH; RV 
COMMON/CLAS/CLASUM < 6 ) 

REAL M(8) 

DATA M(1)«M(2)< M<3). M(4)« M<9)» M(6>* W(7>* M<8) 

«/ 0. Oi 18. 75. 31. 25. 43. 75. 56. 25. 68. 75. 81. 25. 93. 75/ 

C 

C UNPACK ICW TO GET % AREA CODES FOR OPEN WATER < 1 1 ). BARE S0IL<I2). 
C URBANU3), MIXED SOIL AND VEGETATION< 14 > . AND VEGETATION < 15 ) 
I5=ICW/4096 
ICW-ICW-I5M096 
I4=ICW/512 
ICW-ICW-I4*512 
I3==ICW/64 
ICW=ICW'-I3#64 
I2--ICW/8 
I1-ICW-I2*B 
C 

C ZERO BRIGHTNESS TEMPS FROM LAST CALLCIF ANY) 

BTVg=0. 0 
BTHW=0. 0 
BTVB=0. 0 
BTHB=0. 0 
BTVU=0. 0 
BTHU=0. 0 
BTVM=0, 0 
BTHM*0, o 
BTVV=0. 0 
BTHV=0, 0 
BTVF=0. 0 
BTHF=0, 0 
C 

C CALL SUBROUTINES IF CODE WORDS II. . . 15 ARE NOT =0 
IF<I1. NE, 0)CALL WATER (NB. TP, BTVW, BTHW) 

IF< 12. NE. 0)CALL BARE(NB, TP, SM, ROU, BTVB. BTHB) 

IF< 13. NE. 0)CALL URBANCTP. BTVU. BTHU) 

IF( 14, NE, 0)CALL MIX(NB. TP. SM. ROU, BTVM. BTHM) 

IF< 15. NE. 0)CALL VEG<NB, TP, SM. BTVV, BTHV) 

C 

C DETERMINE IF SUBROUTINE FOREST NEEDS TO BE CALLED 
SUM=W(I1 + 1)+W(I2+1)-^W< 13+1 )>W(I4-*-l)+W( 15+1) 

WFbR=100. -SUM 
IF<WFOR. LT, 0. >WF0R=0. 0 
C SUM 7. OF EACH CLASS IN EACH UNIT CELL 
CLASUM( 1 )*W< 1 1+1 >+CLASUM( 1 ) 

CLASUM ( 2 ) «W ( I 2+ 1 } +CLASUM ( 2 ) 

CLASUM ( 3 > *W ( 1 3+ 1 ) +CLASUM ( 3 ) 

CLASUM ( 4 ) »W ( 1 4+ 1 ) +CLASUM < 4 ) 

CLASUM < 5 ) -W ( 1 5+ 1 ) +CLASUM ( 5 ) 
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CLASUH ( 6 ) *WFOR+CLASU« < 6 ) 

IFIWFOR. GT. 5. )CALL FOREST (TP* DTVF* DTHF) 

C 

C calculate composite BRIGHTNESS TEMPERATURES 

SUMBTV*W(1 1 + 1 >*DTVW"»-W( 12+1 )*BTVB+W< 13+1 )*BTVU+W< 14+1 )#BTVM 
«+W( 15+1 )*BTVV+WFOR*DTVP 

SUHBTH»Wm + l >#BTHW+W(12+1 )*BTHB+W< 13+1 )+BTHU+W< 14+1 )#BTHM 
«+W < 1 5+ 1 > #BTHV+WFOR#BTHF 
BTU50=SUMBTV/100. 

BTH50*SUMDTH/100. 

RETURN 

END 
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SUBROUTINE MATER <NB, IF, DTVM, BTHM) 

C CORRECT MATER TEMP <TMC< DEC C) — TMK(DEG K)> 
TMC=(TP-’25. )#0. 254-25. 

TMK-TMC+273. 15 

ROUTE ACCDNQ TO DAND#^ ND 
CO TO < 100, 200, 300), NB 

CALCULATION FOR L-BAND 
100 EH50=0. 256-4TMC*0. 000467 
EU50*0, 5054-TMC*0. 000767 
SKYT*6. 0 
GO TO 400 

C-BAND 

200 EH50=0. 265 
EV50=0. 522 
SKYT=8. 0 
GO TO 400 

X-BAND 

300 EH50=0. 28B-TMC*0, 0003 
EV50=0. 557-TMC*0, 0005 
SKYT=10. 0 
C 

C COMPLETE CALCULATIONS 

400 BTVM=EV50^»TMK+< 1. -EV50)*SKYT 
BTHM*EH50*TMK+( 1, -EH50)*SKYT 
RETURN 
END 
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SUBROUTINE URBAN ( TP . BTVU, BTHU ) 
C CALCULATION FOR ALL BANDS 
TSK=TP-f273. 15 
BTVU*TSK#0. 96 
BTHU»TSK*0. 86 
RETURN 
END 
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SUBROUTINE WATER (ND, TPi BTVW* BTHW) 

C CORRECT WATER TEMP(TWC<PEG C ) --TWI^ < DEO K>) 
TWC=(TP-25. )* 0 , 23+25. 

TWK-TWC+273. 13 

ROUTE ACCDNG TO DAND#i ND 
00 TOUOO, 200, 300), NB 

CALCULATION FOR L~BAND 
100 EH30*0. 256+TWCtO. 000467 
EU30«C. 50S+TWC#0. 000767 
SKYT»6. 0 
GO TO 400 

C-BAND 

200 EH50*0. 265 
EV50»0. 522 
SKYT*8. 0 
GO TO 400 

X“BAND 

300 EH50»0. 288“-TWC#0. 0003 
EU50«0. 557-TWC+O. 0005 
SKYT-IO. 0 

COMPLETE CALCULATIONS 
400 eTVW~EV50*TWK+(l. •“EU50)+SKYT 
DTHW«EH50#TWK+ ( I . -EH50 > ^»SKYT 
RETURN 
END 
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SUBROUTINE BARE(NB< TP» SNi ROU, BTVB» BTHB) 
COMMON/ROUGH/RH, RV 

ROUTE ftCCDNG TO BAND#, NB 
GO TOC 100, 200, 300), NB 

L-BAND ALGORITHM 

100 TGK=250. 15-0, 26*Sh+TP 
IFCSM. GT. 38. >TGK=TP+240, 15 
IFCSM. GT. 12, )G0 TO 101 
EH50=0. 9 - 0 . 00917*SM 
EV50=0. 98-0. 0025*SM 

GO TO 500 

101 EH50*0. 96-0. 0139#SM 
EV50=1. 047-0. OOSOe#SM 
GO TO 500 

C-BAND 

200 TGK=260. 15-0. 53»SM+TP 
IFCSM. GT. 38. )TGK=TP+240, 15 
IFCSri. GT. 12. )G0 TO 201 
EH50=0, 86-0. 00833*SM 
EV50=0. 97-0, 0025*SM 

00 TO 500 

201 EH50=0. 92-0. 0135<»-SM 
EVH50=1. 04-0. 00846*SM 
00 TO 500 

X-BAND 

300 TGK=273. 1 5+TP-O, B7-aSM 
IFCSM. GT. 38. )TGK=TP+240. 15 
IFCSM, GT. 12. )G0 TO 301 
EH50=0. 91-0, 00917-K-SM 
EU50=0, 99-0. 0025*BM 
GO TO 500 

301 EH50=0. 96-0, 0135*SM 
EV50-1. 05-0. 0077-«-SM 

COMPLETE CALCULATIONS 
500 RH=1.-EH50 
RV=1. -'EV50 

RFAC=EXPC-R0U*0. 4132) 

RH=RH*RFAC 

rv=rv*rfac 

EH50=1. -RH 
EV50=1.-RV 
BTHB=EH50*TGK 
BTVB=EV50*TGK 
RETURN 
END 
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SUBROUTINE MIX(NB, TPi SM, ROUi BTVM* BTHM) 

COriflON/ROUGH/RH* RV 

CALL BARE(NBi TP» Sf1» ROU# BTVB. DTHB ) 

CALL VEG(NB* TPi SM» BTW. BTHV) 
BTVt1*(BTVV+BTVB)»0. 5 
BTHM»<BTHV-*-BTHB)#0. 5 
RETURN 
END 
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SUBROUT I NE VEG ( NB , TP , SM. BTW, BTHV ) 

COMMON/ROUGH/RH» RV 

ADJUST CANOPY TEMP* TVC 
TVC»< TP-25. )*0. 25+25. 

TVK*TVC+273. 15 

CALL BARE FOR BASIC EMISSION DATA 

CALL BARE<NB» TVC, SM, 0. 0, BTVB, BTHB ) 

COMPUTE VEGETATION CORRECTION FACTOR ACCDNG TO SOIL MOISTURE 
VFAC=0. 8-0. 00395+SM 

APPLY VEGETATION CORRECTION FACTOR 
XRH=RH*VFAC 
XRV=RV#VFAC 
S0LTMP=BTVB/<1. -RV) 

BTVV=< 1. -XRV)*SOLTMP 
BTHV=( 1. -XRH)+SOLTMP 

RETURN IF L-BAND 

IF(NB. EQ, 1 >RETURN 

X-BAND CALCULATION(CONSTANT EMISSIVITY FOR CANOPY) 

TV=TVK*0. 95 
TH-TVK*0. 92 

BRANCH IF C-BAND 

IF(NB. EQ. 2) GO TO 200 

BTVV=TV 

BTHV=TH 

RETURN 

200 BTVV=(BTVV+TV)-»0. 5 
BTHV=(BTHV+TH)*0. 5 
RETURN 
END 
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FUNCTION BRTMPdA. TBi P) 
BRTMP*?TA#COS ( P ) **2+TB#SI N ( P ) **2 
RETURN 
END 
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SUBROUTINE BCORR(EVSO/ EHSOi ANOLEi BTVi DTH> 


REAL 

FH(81), FV(81 > 






DATA 

FV( 1), 

FV( 2), 

FV< 3), 

FV< 4), 

FV( 5), 

FV< 6) 

e/ 

. S39745, 

. 539915, . 

540428, . 

541283, , 

542481, 

. 544022/ 

DATA 

FV< 7), 

FV< 8), 

FV( 9), 

FV( 10), 

FV( 11 ), 

FV< 12) 

0/ 

. 545907, 

. 548136, . 

550713, . 

553637, , 

556910, 

. 560535/ 

DATA 

FV( 13), 

FV< 14), 

F V ( 1 5 ) ^ 

FVtl6), 

FV< 17), 

FVC 18) 

0/ 

, 564511, 

. 568844, . 

573533, . 

578682, , 

583973, 

. 589768/ 

DATA 

FV<19), 

FV<20), 

FV(21 ), 

FV<22), 

FV(23), 

FV<24) 

0/ 

. 595911, 

. 602424, . 

609311, . 

616574, . 

624216, 

. 632241/ 

DATA 

FV(25>, 

FV<26), 

FV(27), 

FV<28), 

FV(29), 

FV(30) 

0/ 

. 640652, 

. 649453, . 

658646, . 

668234, . 

678221, 

, 688610/ 

DATA 

FV(31), 

FV(32), 

FV(33), 

FV(34), 

FV(35), 

FV<36> 

0/ 

. 699404, 

. 710604, . 

722215, . 

734237, . 

746672, 

. 759522/ 

DATA 

FV(37), 

FV ( 38 ) , 

FV(39), 

FV(40), 

FV(4l ), 

FV ( 42 ) 

0/ 

. 772787, 

, 786467, . 

800562, . 

815070, . 

829987, 

. 845311/ 

DATA 

FV<43), 

FV<44), 

FV(45), 

FV<46), 

FV<47), 

FV<48) 

0/ 

. 861037, 

.877156, . 

893661, . 

910541, . 

927783, 

. 945371/ 

DATA 

FV(49>, 

FV(50), 

FV<51 ), 

FV(52), 

FU<53), 

FV<54) 

0/ 

, 963285, 

. 981504, 1. 

, 1. 

018742, 1, 

037694, 

1. 056812/ 

DATA 

FV<55), 

FV<56), 

FV(57), 

FV(58) 



0/ 

1. 076046. 1. 095339, 1. 

114622, 1. 

133818/ 



DATA 

FU ( 59 ) . 

FV<60), 

FV<61 ), 

FU<62) 

FV(63) 

0/ 

1. 152837, 

1. 171573, 

1. 189907, 

1.207700, 1.224790/ 

DATA 

FV<64), 

FV<65), 

FV<66), 

FV(67) 

, FV(68) 

0/ 

1. 240991, 

1. 256089, 

1. 269835, 

1. 281941 

, 1,292072/ 

DATA 

FV(69), 

FV(70), 

FV(71 ), 

FU<72) 

, FV(73) 

0/ 

1. 299842, 

1. 304801, 

1. 306423, 

1.304097, 1.297109/ 

DATA 

FV<74), 

FV(75), 

FV(76), 

FV<77) 

, FV(78) 

0/ 

1. 284619, 

1. 265646, 

1. 239028, 

1,203400, 1.157136/ 

DATA 

FV(79), 

FV<80), 

FV(81 ) 




0/ 

1. 098311, 

1, 024620, 

0. 933305/ 



DATA 

FH< 1), 

FH( 2), 

FH( 3), 

FH( 4), 

FH< 5), 

FH< 6) 

0/ 

, 539745, 

. 539574, . 

539061, 

538206, . 

537007, 

. 535463/ 

DATA 

FH< 7), 

FH( 8), 

FH( 9), 

FH(IO), 

FH ( 1 1 ) , 

FH< 12) 

0/ 

. 533572, 

. 531333, , 

528742, 

525796, , 

522494, 

. 518829/ 

DATA 

FHU3), 

FH(14), 

FH<15), 

FH(16), 

FH(17>, 

FH< 18) 

0/ 

. 514799, 

. 510399, 

505623, 

500466, , 

494922, 

. 488986/ 

DATA 

FH( 19),. 

FH(20), 

FH(21 ), 

FH(22), 

FH(23), 

FH(24) 

0/ 

. 482649, 

. 475906, 

468747, 

461165, . 

453150, 

. 444695/ 

DATA 

FH(25), 

FH(26), 

FH(27), 

FH ( 28 ) , 

FH ( 29 ) , 

FH<30) 

0 / 

. 435788, 

, 426421, 

416581, 

406258, . 

395440, 

. 384113/ 

DATA 

FH(31). 

FH(32), 

FH(33), 

FH(34), 

FH<35), 

FH(36) 

0 / 

. 372266, 

. 359884, 

346953, 

333459, . 

319385, 

. 304715/ 

DATA 

FH(37), 

FH<38), 

FH(39>, 

FH(40), 

FH<41 ), 

FH(42) 

0 / 

. 289434, 

. 273522, 

256961. 

239734, . 

221819, 

. 203198/ 

DATA 

FH < 43 ) , 

FH<44), 

FH(45>, 

FH(46), 

FH(47), 

FH<48> 

0 / 

. 183846. 

. 163745, 

142870, 

121198. . 

098705, 

, 075365/ 

DATA 

FH(49), 

FH(50), 

FH(5l), 

FH(52), 

FH(53). 

FH<54) 

0 / 

. 051152, 

. 026040, 

000000 , - 

026996. 

054978. 

-. 083974/ 
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SUBROUTINE C ALT (RESL, ALTO) 

COMMON/GAN/BWIDi XPI» IF 
COMMON/OR lEN/ROLL, AZIM* ALT* A. B 

XJ=TAN< (ROLL+BWID/2. >/S7. 3)-TAN< (ROLL-BWID/2. >/57. 3) 
XM*TAN< <BWID/2. )/57. 3) /C0S<R0LL/»7. 3) 
ALT0»RESL/SQRT<2. *XJ*XM) 

RETURN 

END 


